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ABSTRACT 


Laminated fiber reinforced composite materials such as Graphite/ 
Epoxy are generally designed using elastic considerations. Although 
graphite fibers are essentially elastic, the epoxy matrix behaves 
in a viscoelastic manner. The resulting Graphite/Epoxy composite 
material exhibits creep and delayed failures. Time dependent processes 
which are quite slow at room temperature are accelerated by higher 
temperatures and other factors. Assuming the applicability of the 
Time Temperature Superposition Principle (TTSP) concept, short term 
experimental creep compliance and creep rupture data should be useful 
in predicting the long term behavior of laminates at lower temperatures. 
Such an accelerated characterization procedure should have an impact 
on the design of laminated composite structures where combinations of 
temperature, moisture content, applied stress level, and duration of 
load application may necessitate the use of a time dependent analysis. 

An incremental numerical procedure based on lamination theory 
is developed to predict creep and creep rupture of general laminates. 
Existing unidirectional creep compliance and delayed failure data 
is used to develop analytical models for lamina response. The 
compliance model is based on a procedure proposed by Findley which 
incorporates the power law for creep into a nonlinear constitutive 
relationship. The matrix octahedral shear stress is assumed to control 
the stress interaction effect. A modified superposition principle is 
used to account for the varying stress level effect on the creep strain. 
The lamina failure model is based on a modification of the Tsai -Hi 11 



theory which includes the time dependent creep rupture strength. A 
linear cumulative damage law is used to monitor the remaining lifetime 
in each ply. 

Creep compliance and delayed failure data is presented for 
several general laminates along with the numerical predictions. 

Typical failure zone pictures are also given. The compliance predic- 
tions for matrix dominated laminates indicate reasonable agreement 
with the experimental data at various stress levels. Predictions for 
fiber dominated laminates are erroneously bounded by lamination 
theory assumptions. Failure predictions are of the right magnitude 
but are not in exact agreement. Reasons for these discrepancies are 
presented, along with recommendations for improving the models and 
the numerical procedure. 
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INTRODUCTION 

The use of fiber reinforced materials is a concept that dates 
back at least to the use of straw in sundried Egyptian bricks. In the 
past two decades, however, there has been intense new interest in the 
use of relatively strong, stiff fibers as reinforcement in an other- 
wise weak and compliant matrix. Glass, boron, graphite, kevlar, and 
other fibers have been used in either chopped or continuous filament 
form with a variety of matrix materials including epoxy, polyester and 
aluminum. 

Chopped fiber composites have been injection molded to form 
panels and more intricate shapes such as automobile grills. Continuous 
filaments impregnated with resin have been wound around mandrels to 
produce lightweight pressure vessels, missile cases, and struts for 
spacecraft, as well as more domestic items such as golf clubs, fishing 
rods, and bicycle frames. Continuous filaments may be arranged in uni- 
directional plies or woven into a coarse cloth, each of which are then 
impregnated with resin. These prepreg laminae may be stacked at 
various fiber angles and thermally cured in autoclaves to produce 
stiff, lightweight panels, spars, and fairings. Fabrication of other 
structural components by other techniques is also possible. 

The term "advanced composites" has been applied to such continuous 
filament systems as boron/epoxy and graphite/epoxy (Gr/Ep) to dis- 
tinguish them as much stronger and stiffer materials than other 
composite material systems. Advanced composites provide the strength 
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and stiffness of structural metals at a fraction of the weight. 

Primarily because of cost and performance among all continuous fila- 
ment composites, graphite/epoxy currently finds the most wide-spread 
use and is the material of primary interest in the present study. 

While the current applications of graphite/epoxy are primarily 
limited to high performance military aircraft and spacecraft, and a few 
consumer products in the area of sports equipment, there is much 
interest in introducing this material into other areas. Possible appli- 
cations to the transportation industry support substantial current 
interest. The great potential for composite materials is derived from 
its reduced weight, improved fatigue resistance, greater design flexi- 
bility for tailoring material properties to meet design requirements, 
reduced manufacturing costs and fabrication scrap, and improved dimen- 
sional stability due to lower thermal expansion. There are, however, 
many unknowns still concerning the use of this "new" material system. 

The current literature abounds in work in the area of characterization, 
analysis, and design of composite materials. 

Our interest lies in studying time dependent viscoelastic stress- 
strain response of polymer based composite materials and developing 
techniques to predict the long term creep and creep rupture properties 
based on short term testing. Such accelerated characterization T . 

procedures are of obvious practical value to the designer, who cannot 

Vi 

afford to wait for the results of a 10 year test in designing a product 
for a similar intended service life. While graphite fibers have been 
shown to be essentially elastic, the epoxy matrix, as with most 
polymers, exhibits significant viscoelastic response [38,82]. Because 
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the response of a general laminate is governed by the matrix properties, 
as well as the fiber properties, it is important that the laminate be 
considered as a yiscoelastic material. 

For many applications, composite structures can still be 
designed with linear elastic analysis. However, there are applications 
where environmental effects and duration of applied loads require that 
the viscoelastic aspects of the material response be taken into 
account to insure long term structural integrity [45a, b,c]. The 
identification and analysis of these situations provided the impetus 
for the current investigation. 

Previous Efforts 

The present work is a continuation of a. col laborative effort with 
the Materials Science and Applications Office of NASA - Ames Research 
Center and the ESM Department of Virginia Polytechnic Institute and 
State University. The thrust of the project has been to develop tech- 
niques which can give long range strength predictions for general 
laminates. While most of the work deals with graphite/epoxy, it is 
expected that the procedures developed will be applicable to laminated 
composites made with other material systems as well. 

An accelerated characterization was proposed by Brinson [9] of 
VPI & SU to predict the long term response of general laminated 
composite materials based on a minimal amount of material testing. The 
procedure was based primarily upon the time-temperature superposition 
principle (TTSP) which utilizes short term data to predict long term 
results. The work at VPI & SU has been to pursue the development of 



this method to verify assumptions made, obtain data to use with the 
technique, and to determine and correct any problem areas with the 
procedure. A great deal of data has been collected for the graphite/ 
epoxy material system, and substantial progress has been made in time 
dependent characterization. 

The VPI & SU investigations have been directed at visco- 
elastic behavior and creep ruptures, while the NASA - Ames counterpart 
has studied the effects of environment and fatigue on Gr/Ep materials. 
See [63,74,75]. 

Because the VPI & SU work spans several years of research and 
several different batches of material, much of the data previously 
obtained was not directly applicable to the current material. Because 
of variations from batch to batch, and because the properties within a 
given batch have been shown to be highly dependent on the thermal con- 
ditioning [38], there is not as much applicable data available as might 
be expected from such extensive testing. The author's contention is 
that many of the problems encountered in characterizing composite 
behavior are due in part to the variability of the properties from one 
batch to the next--even when made of the same materials and by the 
same manufacturer. The current work is no exception and this has posed 
a great deal of difficulty in interpreting the existing data. 

The Accelerated Characterization Procedure 

The accelerated characterization procedure proposed by Brinson 
[9] is summarized in Fig. 1.1. To characterize a new orthotropic, 
viscoelastic material system, a limited number of tests would be 




Fig. 1.1 Flow chart of the proposed procedures for laminate accelerated 
characterization and failure prediction. 
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conducted with the material to determine ramp loaded static moduli 
and strengths (A). Creep tests would also be conducted to determine 
an Eg master curve and shift function as a function of temperature 
(D). Transformation equations could be used to transform the moduli 
in the material principal coordinate system to any arbitrary fiber 
angle (B,F). It has been found [82] that the shift function is es- 
sentially independent of the fiber angle (E). A time independent 
failure theory can be used to predict static ramp loaded strength of a 
lamina under an arbitrary stress state (C). Based on the assumption 
that the strength master curves have the same shape as the moduli 
master curves, strength master curves may be generated for arbitrary 
stress states and temperatures (G). An incremental lamination theory 
approach would be developed to incorporate the measured lamina proper- 
ties into an analysis procedure capable of predicting the time dependent 
behavior of a general laminate at an arbitrary temperature and subject 
to a given stress state (H). Thus with only a minimum of testing of 
a material system, it is expected that long range predictions of 
strength and compliance for general laminates could be made. Finally, 
long term testing should be performed to verify the validity of the 
procedure (I). 

Many of the ideas incorporated in Fig. 1.1 have already been * 

verified and a great deal of data has been gathered for the graphite/ 

t. 

epoxy unidirectional material. The original work for the accelerated 
characterization procedure can be found with supporting data in 
Brinson, Morris, and Yeow [9] and Yeow, Morris, and Brinson [84]. 

Yeow and Brinson [83] have reported on a comparison of shear 
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characterization methods. Morris, Yeow, and Brinson [55] have 
reported the viscoelastic behavior of the principal compliance 
matrix of the Gr/Ep lamina. Griffith, Morris, and Brinson have in- 
vestigated the nonlinear aspects of the creep compliance [39] and 
creep rupture of unidirectional laminates [40]. 


Outline of Current Efforts 

Much of the characterization of lamina properties has been com- 
pleted in prior efforts as discussed above. A primary focus of the 
current work was an attempt to integrate this unidirectional informa- 
tion into an analysis procedure for a general laminate. A numerical 
method was developed to predict the compliance of a general laminate 
based on the nonlinear compliances of a single lamina. The Tsai -Hi 11 
failure theory was modified for time-dependent strengths and used to 
predict delayed ply failures in the laminate. Experimental compliance 
data for several laminates was taken to investigate the nonlinear 


characteristics of laminates and to 


check the validity of the 


numerical procedure. Delayed failures were obtained for each of a 
variety of laminates tested. Compliance and failure predictions were 
compared with the experimental data. 

Chapter 2 discusses some background concepts and assumptions 
used in the current analysis. Chapter 3 details the development of 
a compliance model as used in the numerical procedure. A discussion 
of the failure model is found in Chapter 4. Chapter 5 provides a 
development of the numerical procedure. Presentation of the experi- 
mental technique is given in Chapter 6, along with details of the 



material system. Chapter 7 expresses the results and comparisons 
of the experimental and numerical investigations. Conclusions and 
recommendations are found in Chapter 8. 



P U ti4- vt O 

onapucr L 


BACKGROUND INFORMATION 

Comments on Terminology 

Prior to discussing the main features of the present endeavor, 
it is worthwhile to clarify the terminology used herein. The term 
laminate refers to the bonded assemblage of several single plies or 
laminae into a panel. Laminate properties refer to the properties of 
the assemblage, while lamina properties refer to the properties of a 
single ply. For practical reasons, a single ply would be very diffi- 
cult to test. Thus lamina properties are determined from testing 
unidirectional laminates, composed, in our case, of sixteen .0052" 
thick plies. Lamina properties are assumed to be equivalent to those 
obtained from a unidirectional laminate. 

Because the fibers are much stiffer than the matrix, the fiber 
properties tend to control the response of a lamina in the fiber direc- 
tion. Thus the compliance in the fiber direction is said to be a 
fiber dominated property. Also, in a uniaxial test along the fiber 
direction, the transverse strain is closely tied to the axial strain. 
Thus the Poisson's ratio effect, or S-^ term of the compliance matrix 
is also considered fiber dominated unless there has been significant 
degradation of the matrix. On the other hand, the compliances in 
shear and transverse to the fiber direction are referred to as 
matrix dominated properties because they are closely tied to the 
matrix response. 
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The coordinate system convention used for fiber reinforced 
lamina is well standardized and is illustrated in Fig. 2.1. The x-y 
coordinates are referred to as global coordinates. The 1-2 coordinates 
are referred to as the local coordinates or the principal directions 
of the lamina and the 1 direction is parallel to the fiber orientation. 

Balanced laminates are those which for every ply at an angle 
0 to some reference axis, there is also an identical ply at an angle 
-e. These reference axes are the principal axes of a balanced 
laminate. Symmetric laminates are those in which the laminae orienta- 
tions form a mirror-image about the laminate midplane. 


Orthotropic Constitutive Relations 

The primary application of laminated composites is as flat or 
shallow panels loaded by in-plane loads in a state of plane stress. 
As such, the linear elastic constitutive relation for an anisotropic 
material may be simplified from the most general expression. 


e ij S ijkt a k£ 


i , j,k,Jt = 1 ,2,3 


( 2 . 1 ) 


£.. = strain tensor 
S i j k £ = ^ term compliance tensor 
= stress tensor 

to the reduced form applicable to orthotropic lamina under plane 
stress situations = t 13 = 0) [48], 
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( 2 . 2 ) 






12 


or 

{e} = [S]{a> 

where 

e-j , £2 = in-plane normal strains 

y 12 = 2e ]2 = ^ n_ P^ ane engineering shear strain 

S = reduced compliance matrix 

a-| ,©2 = in-plane normal stresses 

t 12 = in-plane shear stress 

The components of the compliance matrix may be expressed in terms 
of the engineering material constants as, 

S 11 = 1/E 11 
S 22 = 1/E 22 

S 1 2 = S 21 = " v 12 /E 11 = “ V 21 /E 22 
S 66 = 1/G 12 

The strains may be transformed between local and global co- 
ordinates by, 

{e> 12 = [T 2 ]{e> xy or { e ) xy = [T 2 ] _1 {e} 12 (2.3) 

and the stresses by 

{a} 12 = [T r ]{a} xy or {a> xy = [T^'^a}^ (2.4) 

where the transformation matrices are given by. 
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The compliance matrix developed in the material principal 
coordinates can also be transformed to the global coordinates, 

fe>xy= C^]- 1 U} !2 = [T 2 r 1 [S][T 1 H° xy > 

or [S] = [T-.r’mtr,] (2.5) 

Alternatively, the constitutive relations may be expressed in 
terms of a reduced stiffness matrix 

{a) 12 = [Q]{ £ } ]2 (2.6) 

where [Q] = [S]~^ 

Similarly, the stiffness matrix can be transformed to the global 
coordinate system by 

[Q] ■ [T i r 1 [Q][T 2 ] (2.7) 

Experimentally, and are obtained from uniaxial tests 
(normally tension tests) on specimens cut parallel and perpendicular 
to the fiber direction of a unidirectional laminate, respectively. 

S-ji and S 22 are determined from axially mounted strain gages (or other 
deformation measuring devices). S-j 2 (= ) may be determined from a 

transversely mounted gage on either specimen. Sgg has been determined 
by a variety of techniques including rail shear, picture frame speci- 
mens, and off-axis tensile specimens. Chamis and Sinclair [16] have 
proposed the use of a 10° off-axis unidirectional specimen to measure 
the shear compliance. Yeow and Brinson [83] made a study of several 
methods for determining shear properties and have concluded that this 
10° specimen is the best configuration for measuring shear properties. 



Use of an electrical strain gage rosette on the off axis specimen 
allows determination of the shear strain. The shear compliance may 
then be computed directly. Alternatively, S^, , as determined from an 
axially mounted strain gage on an off axis specimen, may be used to 
calculate the shear compliance. Expanding the first term of Eq. 2.5 
yields 


§n = S xx = cos^e S-j i + sin^e + cos^e sin^e (2S 
Knowing S-^, one can solve for the shear compliance, S^g, 
S -j -j , S ^ ^ 5 ^12 5 an d ^ • 


12 + S 66^ 


in terms of 


Lamination Theory 

Classical laminated plate theory is an important tool in the 
analysis of laminated composite materials. The basic assumption of 
this theory is that lines normal to the laminate mid-plane remain 
straight and normal to the mid-plane after loading. This implies that 
there are no interlaminar shear deformations or stresses. While this 
is a valid assumption for interior regions of well bonded panels, it 
cannot be physically correct near free edges of the plate where 
interlaminar stresses must occur to maintain equilibrium. However, 
it may be shown that these regions are very localized [48]. The 
assumption of no interlaminar deformations may break down in specimens 
which have undergone large deformations or when failure is imminent. 
This theory is widely used, however, and it was felt that it would be 
adequate for the present analysis. Because all the laminates studied 
herein are symmetric about the mid-plane, only in-plane deformations 
result from in-plane loads and vice versa. Thus, only the in-plane 
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stiffness matrix is developed because out-of-plane deformations and 
loads are not considered. 

To compute the constitutive properties of a laminate, the stiff- 
ness matrices for each contributing ply are transformed into the 
global x-y coordinates and combined to provide the total laminate 
stiffness 

K . 

[A] = Z [Q] K t. (2.8) 

k=l K 

where [A] is the laminate stiffness matrix 

K is the number of plies 

- k th 

[Q] is the laminate stiffness of the k ply in 

global coordinates 

■f" h 

t^ is the thickness of the k un ply 
The elastic laminate strains {e} e and the force resultants { N > 
may then be related as, 

{N> xy = [A](e} xy (2.8) 

or 

tcl x, - rcr 1(N) xy 

J_ J- 

To calculate the ply stresses in the k Ln ply, 

= [Q] k {{e}* 2 - {e}^} = [Q] k {e}® 2 (2.9) 

where {e}| 2 - total laminate strain in 1-2 coordinate system of ply k 
{£ }[ 2 - residual laminate strain in 1-2 system 
where the residual strains are any non-elastic strains such as thermal 
strains, hygroscopic strains, creep strains, etc. 
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Linear Viscoelasticity 

For linear elastic materials, the constitutive equation is given 
by. 


e ij S ijk£ a k£ i,j,k,£ - 1,2,3 

For linear viscoelastic materials under creep loading, the compliance 
can be generalized to a function of time 

= s i<k (*! V < 2 - 10 > 

where 


a k£ (t) = 5 k£ H(t) (2jl) 

and H(t) is the Heaviside function. 

For more general loading states, one may express: 

E ij (t) = s ijkt (t > °k» 0 H(t > + s ijkn (t - H> CT k^ H(t - V 
+ ... 


which may be generalized to the following Duhammel integral: 


■uw ■ 


Sijkitt - 0-f— dx 


( 2 . 12 ) 


This expression is often referred to as the Boltzman Superposition 
Principle and is a consequence of and is only valid for a linear 
material [20]. 

As with the linear elastic case, the viscoelastic compliance 
tensor is symmetric. Schapery [66] has verified this analytically as 
long as each of the constituent phases is symmetric. Morris, Yeow, and 
Brinson [55] have shown that this is borne out experimentally. 


t 


t- 
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For plane stress analysis of an orthotropic material, the compliance 
matrix may be reduced to four independent functions of time [44] 


V,(t) 


S n (t) S 12 (t) 0 


r \ 

°1 

e 2 ( t ) 

■ r 

S-j 2 (t ) S 2 2 ( t ) 0 

. 

s 2 

J i £ ( t ) , 


0 0 Sgg (t)_ 


- T 1 2- 


For our project, all experimental compliance data were obtained from 
uniaxial tension tests. Compliance properties were assumed to be the 
same in compression as in tension. For a linear viscoelastic material, 
Eqn. 2.13 applies to any general plane stress state, although no 
biaxial tests were run for verification. 

An analogous development may be used for the relaxation modulus 

- c iji« (t) ^ (2 - 14 > 

where 

£ k£ (t) = £ k£ H(t) (2J5) 

The creep compliance has been used throughout the present analysis 
because of the difficulty in obtaining a pure relaxation test for the 
current material system. 

Time Shift Superposition Principles 

The underlying premise for an accelerated characterization of a 
material system is that one can in some way use short term experimental 
data to predict long term material response. The Time Temperature 
Superposition Principle (TTSP) has found wide use in polymeric studies 
since its introduction by Leaderman in 1943. The basic idea is that 
compliance curves at different temperatures are of the same basic 
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shape, but only shifted in time. Thus by taking short term compliance 
data at several temperatures and then shifting these curves hori- 
zontally in log time (some vertical shift may also be necessary, see 
Griffith [38]) one can obtain a smooth curve approximating the 
compliance over many decades of time. 

The response of a single Kelvin element (a spring and dashpot 
in parallel) is 

D(t) = E(1 - e~ t/x ) (2.16) 

x = n/E 

where D is the creep compliance, E is the modulus of the spring, x is 
the retardation time, and n is the dashpot coefficient. If N Kelvin 
elements are connected in series, the overall creep compliance is 
given by the Prony series, 

N ., 

D(t) = z E . (1 - e" t/T i) . (2.17) 

i=l 

Any monotonically increasing creep compliance function of a linear 
material may be approximated by a generalized Kelvin model composed of 
many individual Kelvin elements connected in series. As the number of 
Kelvin elements becomes infinite, the fit becomes exact. A finite 
number of elements would result in a discrete distribution of the 
retardation times whereas an infinite number of elements would yield a 
continuous retardation spectrum where the creep compliance may be 
given by. 
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/*°° . i 

Otto.*) - W + L(T o’ £n t)[1 ' e “ ] d t 

+ t/n 0 (T 0 ) (2.18) 

where D is the creep compliance 

T q is a reference temperature 
D q is the initial compliance due to a free spring 
L is the retardation spectrum 
x is the retardation time 

and t/n represents the flow of a free dashpot, which was assumed 
to be zero for the current material 

For Thermorheologically Simple Materials (TSM), the compliance 
at other temperatures is represented by, 

*oo 

D(T,t) =.D 0 (T 0 ) + L(T 0 ,£n x)[l - e" ?/x ]d £n t (2.19) 

where 5 is the reduced time given by, 

5 = t/a T (2.20) 

and a^ is the temperature shift factor. When the retardation spectrum 
and compliance are plotted vs a log time scale, the effect of ay is 
merely to shift these curves to the right or left in time according 
to, 

log 5 = log t - log a y (2.21 ) 

Unfortunately, most engineering materials do not fit into the 
TSM description, and are classified as thermorheologically complex 
materials (TCM). For this case, there will be a vertical shift in the 



retardation spectrum and compliance, as well as a horizontal shift. 
The compliance at a temperature T is given by. 


D(Tjt) = D q (T) + 


L(T,£n t)[1 - e ^ T ]d £n x 


( 2 . 22 ) 


For many materials and temperature ranges, however, the temperature 
dependence of D Q and L tends to be fairly small. The horizontal 
shift for various temperatures remains the fundamental concept. A 
more detailed discussion of these concepts may be found in Ferry [28]. 
To use the TTSP for either TSM or TCM, compliance data is taken for 
a number of different temperatures. The duration of these tests is 
normally quite short because of practical considerations. This short 
term data is then shifted to form a smooth and continuous "master 
curve" which is assumed to be valid over many decades of time at an 
arbitrary reference temperature. Various techniques have been used to 
determine the appropriate amount of horizontal and vertical shift. A 
good discussion of these aspects is found in Griffith [38]. To obtain 
the compliance at other temperatures, the master curve is shifted to 
coincide with the short term data at that temperature. Rosen [62] 
illustrates the application of the TTSP technique for the relaxation 
modulus of a TSM polymer, as is reproduced in Fig. 2.2. 

The point of primary interest to the current paper is simply that 
these techniques have been successfully used to shift compliance data 
obtained at one temperature to predict the compliance at another 
temperature. Shifting procedures are rigorously justifiable only above 
the glass transition temperature, T , although they have been 




TIME (hr) 


Fig. 2.2 TTSP formation process as given by Rosen [62]. 
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successfully employed below the as well. For small stresses and 
strains in the linear range, Yeow [82] found the TTSP to be applicable 
to the current material system. 

Furthermore, others have proposed that delayed yield or failure 
master curves may be constructed in an analogous manner from short term 
data at various temperatures. Lohr [52] constructed yield stress master 
curves for constant strain rate testing of several polymers. Similar 
yield master curves for creep loading were also presented by Lohr, 
Wilson, and Hamaker [53]. No long term verification of their results 
was given. Nonetheless, there does seem to be some justification for 
this technique in the Tobol sky-Eyring Reaction Rate Equation. 

Thus, although most of the current experimental work was con- 
ducted at 320°F, it is possible to utilize such data in predicting 
longer term response at lower temperatures. Once an understanding and 
predictive capability have been established at a given temperature, 
testing at other temperatures can be used to extrapolate this infor- 
mation to long term behavior. 

In addition to temperature, there are also other accelerating 
factors such as moisture content, stress level, cyclic loading, and 
the absorption of jet fuel in wet wing designs. Similar superposition 
principles have been proposed for these factors individually and in 
combinations . 

The effect of moisture on composite laminates is usually pro- 
nounced but is not being studied directly in the current work. Nearly 
constant moisture levels were maintained in the test specimens to 
minimize any influence on measurements. The acceleration due to stress 
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level is of current interest because epoxy behaves nonlinear! y at 
moderate and high stress levels. While moisture effects and cyclic 
loads can be minimized for laboratory experimentation, the stress level 
nonlinearities cannot be avoided. Several nonlinear approaches have 
been proposed by a number of investigators as reviewed by Griffith 
[ 38 ]. 


Unfortunately both the Boltzman and TTSP techniques are referred 
to as superposition principles. Clearly, from all outward appearances, 
the two techniques are unrelated. The first deals with the strain 
resulting from a variable load history and is analogous to the super- 
position principles employed in linear elastic analysis. The latter 
implies that compliance data may be shifted in time and superimposed on 
similar data taken at a different temperature. The former is limited 
to linear viscoelastic materials, whereas the latter applies to a much 
more general response. Conceptually, the two types of superposition 
principles are completely different. It is interesting to note, how- 
ever, that the Schapery procedure described in a following section in- 
corporates elements of both techniques in a nonlinear expression for 
strain due to a general load history. Generalization of this concept 
to include temperature could provide a unified approach to account for 
the two different problems addressed by the individual techniques. 

Nonlinear Viscoelasticity 

A number of techniques have been used to account for nonlinear 
viscoelastic behavior. One such approach is the Time Stress Super- 
position Principle (TSSP) given by the basic equation 
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D(a,T) = D 0 (a) + t> a AD('£) (2.23) 

where 

b is the vertical shift factor 

a 

£ is the reduced time given by £ = t/a 

and a is the horizontal shift factor due to the stress level. 

a 

The similarities between the TTSP and TSSP are apparent. The 
implication is that compliance data at various stresses rather than 
various temperatures may be shifted in log time to predict long term 
compliance based on short term testing. Darlington and Turner [25] 
note that while the TSSP rests on a less rigorous development, it has 
been used with some success. Griffith [38] has done considerable work 
in determining the appropriate horizontal and vertical shift func- 
tions for a combined Time Temperature Stress Superposition Principle 
(TTSSP) application to the current graphite/epoxy material system. 
Griffith's results, however, were not readily adaptable for implemen- 
tation into a numerical scheme. 

The Green-Rivlin Theory, or multiple integral approach has also 
been used to model nonlinear viscoelastic materials. For the one 
dimensional case, the creep strain due to a constant stress is 
assumed to be a polynomial in stress [20]: 

e(t) = D 1 (t)a 0 + D 3 (t,t,t)a 0 3 + D 5 (t,t,t,t,t)a Q 5 + ... (2.24) 

Even powers of a Q are omitted to avoid negative values of the stored 
energy. For general loading, the response is expressed in terms of 
multiple convoluted integrals: 
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D 3 (t-T-|»t T 2 ,t-T 3 ) dx 


da(x-j) da(x 2 ) da(x 3 ) 


dx. 


dx. 


dx^ dx 2 dx^ + ... 


(2.25) 


Arridge [4] notes that the series may be truncated after the third 
order term for some materials, although such a simplification may 


not be accurate in general. He points out that few applications of 
the procedure have been made because of the difficulty in using the 
technique and the prohibitive amount of testing required for general 
characterization. 


Schapery Approach to Nonlinear Viscoelasticity 

Another nonlinear approach of interest is that proposed by 
Schapery [68]. His approach is derived from thermodynamic considera- 
tions and has been used successfully by several investigators [54, 
68,14] to predict the behavior of polymers both with and without fiber 
reinforcement. The form of the constitutive equation for uniaxial 
stress is given by. 


e 


(t) = g_D 
' ' r o o 



dg P o 

AD(ip - \p' ) dx 


(2.26) 


where g , g-j , and g 2 are functions of the streys level, D q is the 
instantaneous compliance, &D is the transient compliance, ip and ip' 
are reduced time parameters as given by. 
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4> = 'i'(t) 



<f>‘ - V' (t) 


fT dr 

o a a 


(2.27) 


where a o is the stress dependent time shift factor. As mentioned 
earlier, the basic form is very similar to the Boltzman superposition 
integral. In fact, in the linear range of the material when a is 
small, g Q = g-| = $2 = a cr = anc * the Boltzman integral is regained. 
Furthermore, the reduced time and shift factor concept is also 
employed. It seems reasonable to hypothesize extending this pro- 
cedure to include the features of temperature superposition by 
perhaps letting 


e(t) - g 0 (<r) D Q (T)o + g^ (a) 


ft d(g ? (a)a) 

AD(T 1 ) ^ — - dx 

— co 


(2.28) 


where ij; and V are now reduced times with respect to both stress and 
temperature shift factors. 


- t(t) 




a T 


ip 1 = Ip 1 ( X ) = 


dt ' 


o 


a a T 

CT T 


(2.29) 


While this approach was not pursued in the current study, the develop- 
ment of a unified technique to account for both temperature and 
nonlinear stress effects would be very advantageous. 

The Schapery procedure is very appealing from the standpoint 
that it provides a unified approach to predicting the nonlinear visco- 
elastic response to an arbitrarily varying stress. The difficulties 
arise in the experimental determination of g , g-| , and a^. 

Because the approach is more general, it requires more information to 


1 


•6 
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evaluate the functions of stress. In particular, Schapery uses creep 
and creep recovery data to determine the unknown functions. Creep 
data alone is insufficient to explicitly characterize the functions 
of stress; all one can obtain is the ratio g^/a^ 

In the previous testing program for the current material system, 
only creep data but no creep recovery data was taken. This limited 
data prevented the utilization of the Schapery procedure at the present 
time. Obtaining sufficient data at various stress levels, temperatures, 
and fiber angles would have been beyond the current scope. 

In discussions with colleagues the authors have been led to 
believe that because the Schapery procedure is so general, the deter- 
mination of unique expressions for g Q , g-j , g^ and a^ is virtually im- 
possible. Apparently, an admissible set of expressions obtained from 
creep-creep recovery, may not be valid for another load history such 
as a multiple step loading. If unique expressions cannot be obtained 
experimentally, the whole procedure will be of little practical use. 

A carefully controlled test program could substantiate these conten- 
tions, or validate the technique for the current material system. 

For the materials he investigated, Schapery proposed that the 
transient compliance would be given by a simple power law which is 
not a function of stress. 

AD ( ip ) = m (2.30) 

By forcing the compliance to be independent of stress, the necessity 
for a stress shift factor is created. The concept of a stress shift 
factor is not required in other nonlinear procedures because the 
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compliance is expressed as a function of stress. For creep loading, 
it can be shown that the two approaches are equivalent. 

Findley Approach to Nonlinear Viscoelasticity 

A nonlinear viscoelastic characterization method extensively 
studied by Findley [30-34] was eventually used in the current analysis. 
The basic concept behind the Findley analysis is that for any given 
creep load, the specimen strain is given by, 

e(t) = e Q + m t n (2.31) 

where s , m, and n are material properties. Further, the assumptions 
are made that 

n = constant, independent of stress level 

e Q = e Q sinh a/a £ (2.32) 

m = m‘ sinh a/a^ (2.33) 

I 

where -e , a £ , m 1 , and a m are material constants for any given tempera- 
ture, moisture level, etc. The nonlinear effect of stress is accounted 
for by the hyperbolic sine terms. Apparently, the approach is 
essentially empirical, although there is some basis in the reaction 
rate equation [30,68]. Nonetheless, Findley's technique was found to 
provide an accurate means to express the current experimental results. 

Time Independent Failure Criteria 

Numerous failure criteria have been proposed and used with vary- 
ing degrees of success to predict static strengths of general laminated 
composites. It is expected that an extension of these strength 
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criteria to include time dependent effects can be used to predict 
creep rupture in general laminates. While there are several basic 
approaches to predicting static strengths, the most widely used method 
independently compares the stress (or strain) state in each ply 
against a lamina failure criteria. If any ply has "failed", the 
properties of that ply are reduced to reflect the damage sustained due 
to failure. If there are intact plies remaining, the load may be in- 
creased and the process repeated until total laminate failure occurs. 

Several of these failure criteria are described below: 

Maximum Stress : The maximum stress failure criteria is a simple, 

straightforward approach involving comparison of the ply stresses in 
principal material directions against their respective critical 
values in tension and compression. 


X < a, < X. 
c I t 


V c < < Y t 


Maximum Strain : The maximum strain criteria is similar to the maximum 

stress theory except ply strains in principal material directions are 
compared against their respective failure strain value 


£ 1 < e l < £ 1, 

c t 


c c t 

I Y 1 2 1 
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While these two theories are easy to apply, they do predict 
cusps in the failure stress vs fiber angle which are not borne out by 
experimental data [77]. 

Tsai -Hi 11 : Hill proposed an extension of the von Mises distortional 

energy yield criteria to anisotropic materials: 

(G + H)a^ + (F + H)a£ + (F + G)^ - 2H o-jc^ - 2G a-jO^ 

- 2F + ^ t 23 + ^ t ^3 + ^ t 12 = ^ (2.34) 

For a plane-stress analysis, this may be reduced to 

■\ 2 r 2 r 

a l °l a 2 a 2 t 12 

X “2 Y S 

J A l J l 

There is no distinction for compressive or tensile critical stress 
values, but investigators have used X c when a-j is compressive, X^ for 
a-j tensile, etc. 

Equation (2.35) is generally accepted to be a more accurate 
representation of experimental data than the previous theories. One 
drawback is that this method only predicts the occurrence of failure but 
does not predict the manner in which failure will occur. 

Tsai - Wu : Another quadratic failure theory is the Tsai - Wu Tensor 

Polynomial criteria. 

F^. + F.ja.aj + ... = 1 i ,j = 1 ,2, . . . ,6 (2.36) 

where F^ and F^ are second and fourth order strength tensors 
respectively. For plane-stress, this failure criteria may be expressed 



2 

=1 (2.35) 



as : 
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F l°l + F 2 a 2 + F 6 t 12 + F ll a l + F 22°2 + 2F l2 a l°2 


+ F 66 t 1 2 1 


(2=37) 


Although quite similar to the Tsai -Hi 1 1 approach, this method is more 
general in the sense that it can account for strength differences in 
tension and compression and provides for independent interactions 
between the normal stress components. 

Note, however, that the independence of the cr-^ interaction 
effect does not permit the accurate determination of F-^ from uni- 
axial tests. The inconvenience of running biaxial tests to determine 
F-J 2 renders this greater generality more of a liability than an asset. 
Because of the increased number of parameters, however, this method 
does tend to be slightly more accurate than the Tsai -Hi 11 formulation. 
The improved accuracy does not usually warrant the extra trouble, 
and Tsai -Hi 11 finds wider use. 

Sandhu Analysis : Another approach to failure criteria is that of 

Sandhu [65]: 


f(cr,e) = K. 


itn i 


- £ i 


°i de i 


= 1 i = 1,2,6 


(2.38) 


where 




L' e i u 


°i de i 


(2.39) 


The appeal of Sandhu 's procedure is that account is made for 
material nonlinearities and failures are based on total energy sus- 
tained by the material. This approach is somewhat inconsistent with 
other failure criteria in that it is based on total energy rather 
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than distortional energy but may have some merit. 

Puppo-Evensen : While the previous techniques involve application of 

the particular failure criteria in a plywise fashion, the Puppo- 
Evenson approach [58] uses a failure criteria based on the laminate as 
a whole. Claim is made that the method incorporates interlaminar 
effects and does not require lamination theory or constitutive equa- 
tions. Admittedly, interlaminar effects are neglected with lamination 
theory approaches; however, there is no rigorous correlation with 
actual interlaminar effects in the P-E theory either. The short- 
coming of this theory is that it is valid only for predicting failure 
due to general loading on the one specific laminate being studied. 
Obviously, such a technique may be quite accurate, but of limited 
usefulness to the designer who has the option to vary the layup. Yeow 
[82] has used the P-E criteria in a plywise manner. In addition to 
being cumbersome to apply, this defeats the purpose of the tensorial 
approach required in the original development of the P-E theory [58]. 

Several basic types of static failure criteria have been men- 
tioned. Other techniques exist, but most require large amounts of 
biaxial data or other properties which are difficult to obtain. An 
excellent review of static, orthotropic failure criteria may be 
found in Rowlands [64]. 

Time Dependent Failure Criteria 

A number of time dependent failure criteria have been proposed 
to predict the time to failure of different materials. Most of these 
techniques are valid only for a uniaxial, constant stress state in 
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homogeneous materials. A relationship often credited to Zhurkov 
has been used quite successfully by Zhurkov [85] to predict the time 
to failure of a wide class of materials. This relationship is 
given by, 


t r = t Q exp 



(2.40) 


where 

t r = rupture time 

t Q = a material constant supposedly based on atomic vibrations. 
Zhurkov contends that this term is the same for most 
materials. 

U Q = activation energy (material constant) 
y = a constant 

a = applied uniaxial true stress 
K = Boltzman's constant 
T = absolute temperature 

While the technique is highly acclaimed in the Russian literature, it 
is not as widely accepted among other investigators. 

Slonimski et al [70] have modified the basic Zhurkov equation to. 


t r = t Q exp 




(2.41) 


This form, known as the modified rate equation, has been successfully 
used by Griffith [38] to fit the delayed failure data of 90°, 60°, 
and 45° off axis specimens at 290, 320, 350, and 380°F. Predictive 
capabilities of the procedure have not been verified, however. 
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A variety of other time dependent stress limit failure criteria 
have also been proposed. A discussion of several of these methods 
may be found in Griffith [38], and an excellent overview of a number 
of these techniques is presented by Grounes [43]. 

Of particular interest in the current analysis is an extension 
of the Tsai-Wu tensor polynomial for anisotropic materials to include 
time dependent strengths. Wu and Ruhmann [81] have proposed this 
technique and applied it to unidirectional glass/epoxy composites. 
They envision a Tsai-Wu static failure surface (t Q ) in a-j , T ]2 

space. Other surfaces within F(t Q ) describe the time dependent 
strength for any arbitrary stress state vector and are given by, 


F(t) = 


£ dT + F(t o> 


The integral reflects the decreasing magnitude of the strength vector 
with time. Wu and Ruhman have suggested that this could be an 
exponential decay following Zhurkov. This is a classic paper con- 
taining statistical analysis of data obtained from room temperature 
creep rupture in an air and a hostile benzene environment. 



Chapter 3 

CONSTITUTIVE BEHAVIOR MODEL 

There are a wide variety of approaches that could be used to 
model the constitutive properties of an orthotropic viscoelastic 
material. The criteria used to select the model subsequently developed 
was for the approach to be nonlinear and to provide a good fit for the 
existing unidirectional compliance data. Also, an important considera- 
tion was for the model to be a fairly simple approach which could 
easily be adapted to the numerical scheme developed in Chapter 5. 

There are several difficulties in extending existing theories to the 
case of a variable, biaxial stress state for a nonlinear viscoelastic 
orthotropic material. These problems and the approach eventually used 
will be discussed in this chapter. 

Constant Uniaxial Stress for Isotropic Materials 

In order to analyze the compliance of the current material 
system, a necessary consideration was to understand how to characterize 
the creep compliance for the simplest case--creep of an isotropic 
material under a constant uniaxial stress. Hundreds of studies have 
been conducted for creep of different materials, different condition- 
ing (e.g., aging), different temperatures, and ways to predict the 
response, the temperature effect, and the nonlinear stress effect. 
Nearly all have only dealt with this simplest situation. It is only 
fitting that the study of a variable biaxial stress state in an 
orthotropic material should begin here. 
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Fessler and Hyde [29] have suggested that much of the work in 
the area of predicting creep has been the characterization of the 
following type expressions for the initial component of strain 

e r = o/E + f-j (ct) f 3 (I) (3.1) 

and the creep strain 

£ c = V 0 ) f 2 (t) f 3 ( T) (3-2) 

The assumption of separation of variables seems to be one of con- 
venience rather than physical reasoning. 

Common types of stress dependence are: 

f-j (a) = Aa" 1 (3.3a) 

f-|(a) = A sinh (a/a Q ) (3.3b) 

or f-j(a) ■= A exp (a/a ) (3.3c) 

The hyperbolic sine expression will subsequently be used in the 
current investigation. It should be noted that this form falls between 
the other two expressions, tending towards Aa m for small a, and 
towards A exp (a/a ). for large values of a. 

Expressions for the time dependence are very numerous. The most 
common is the power law: 

f 2 (t) = t n . (3.4) 

Conway [22] discusses a wide variety of other expressions which have 
been used with varying degrees of success for many materials. These 
range from logarithmic forms to polynomials in time: 

f 2 (t) = at 1/2 + bt + c t 3/2 
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to the famous Andrade one-third creep law: 
f 2 (t) - (1 - Bt 1/3 ) e kt - 1 

Fessler and Hyde [29] point out that the temperature dependence 
is almost invariably assumed to be 

f 3 (T) = exp (-U/KT) (3.5) 

where 

U = activation energy 
K = Boltzman's constant 
T = absolute temperature. 

Supposedly, this expression is fundamental to all rate processes. 

The Power Law for Creep 

A power law representation for transient strain is independent of 
stress or temperature dependence assumptions provided the uncoupled 
form of Eqs. 3.1 and 3.2 is appropriate. Conway [22] points out that 
the power law is by far the most widely used form, and is applicable 
to a wide variety of materials. 

Consider the power law in the form 

£ ( t ) = £ q + mt n (3.6) 

e(t) = nmt n- ^ (3.7) 

where £ Q , m, and n are material parameters valid at a certain stress, 
temperature, etc. In Eqn. 3.6, £ Q is referred to as the initial or 
instantaneous component, and the second term represents the transient 
or creep component of the strain. Note that while s Q is often 
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considered to be the initial or instantaneous strain, it is actually 
a curve fitting parameter. As such, £ Q may not necessarily correspond 
to an actual instantaneous strain even if this value can be physically 
measured. In fact, often the assumption of e Q = 0 is made in cases 
where the instantaneous response is small in comparison to the total 
strain. For these cases, the predicted total strain may provide a 
good fit for the data over a considerable time range, but not be valid 
at very short times [54]. 

Several techniques can be used to determine the material 
constants for the power law from experimental data. Remembering that 
a power law plots a straight line on log-log paper, an obvious 
procedure is to use a trial and error approach for determining e . The 
correct value of e would result in the best fit to a straight line on 
log-log paper of the transient component of strain. The slope of the 
line gives the value of the exponent n, and the t = 1 intercept is 
the value of m. This method, although probably the most accurate, is 
also very tedious. 

Another approach is to record the strains e-j , and 
times t-j , t£, and t^, where 



The power law parameters may be easily determined from the following 
equations as found in Boiler [7]: 





* 


n 


(3.9) 
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£ 


0 


2 

£ ] £ 3 ~ £ 2 
£] - 2s 2 + £3 


(3.10) 



This approach, while much simpler, is probably not as accurate as the 
preceding method because the fit is based on only 3 data points rather 
than a larger number as might be used in a graphical trial and error 
approach. Obviously greater care must be exercised in using this 
simpler method. 

Conway [22] discusses another method to determine the power law 
parameters which is based on creating a point to point difference 
table from the experimental data. Values of m and e Q are then calcu- 
lated for each step. The average value of m and e Q are assumed to be 
the best values for these parameters. This technique is worthy of note 
although the procedure was not used in the current analysis. 

Eqns. 3.9, 3.10, and 3.11 were used in the current study to 
avoid the tedious process of plotting data for various guesses of 
£ q until a satisfactory value had been obtained. The time values t-j 
and t^ were chosen to span the time range for the experimental data. 

Perhaps part of the reason the power law has found such wide 
usage is its great versatility to represent a variety of material 
responses. Depending on the value of the exponent n, the power law 
may be used to describe several viscoelastic material types as is 
shown in Fig. 3.1. 



LOG TRANSIENT STRAIN 
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LOG TIME 

Fig. 3.1 Versatility of the power law in predicting a variety of 
material responses. 
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Both e(t) and e(t) approach infinity as t increases without 
bound for n > 1. This region has been labeled as a super fluid in 
Fig. 3.1 because the strain rate at a constant creep load increases 
with time. While not relevant to the present materials, this region 
could be used to characterize fluids with a decreasing viscosity such 
as perhaps an engine oil deteriorating with usage. Such response is 
known as shear thinning or pseudoplastic [51], 

For the case of n = 1 , the response is that of a viscoelastic 
fluid. Specifically, this represents a Maxwell element which may 
consist of a nonlinear spring and/or dashpot. Despite this generality, 
the strain rate is always a constant, e, at any given stress, thereby 
limiting the usefulness of this equation for real fluids. 

The region of 0 < n < 1 accounts for most practical applications 
of the power law. As n approaches a value of unity, the response is 
fluid-like; as n approaches zero, the behavior is solid-like. For 
intermediate values, however, the behavior is neither that of a true 
solid, because the strain increases without bound, nor of a true fluid, 
because the strain rate approaches zero. Actually, this is the 
accommodating feature of the power law because most engineering 
materials are neither true fluids nor true solids, but somewhere in 
between. 

For n = 0, the obvious conclusion is that the response is not 
time dependent. However, as will be discussed later, Eqns. 3.9, 

3.10, and 3.11 can predict singular values of e Q and m when n = 0. 
Singular values of and m may suggest a time dependent response for 
this case but bounded values yield time independent response. 
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While most power law studies require that m and n be positive 
quantities, interestingly the case of n < 0 actually describes a true 
viscoelastic solid. This condition simply implies that the strain 
is a bounded quantity and requires that e Q > e(t) and m < 0. 

One drawback is that the power law does not have a simple 
mechanical analog as does the generalized Kelvin element, for instance. 
The lack of a mechanical analogue poses no problem mathematically but 
one may miss the physically meaningful features of a system of linear 
springs and dashpots. Lockett [51] contends that the power law 
represents a simpler approach than the use of several exponential terms 
associated with a generalized Kelvin element, and is often just as 
accurate. His reasoning is quite understandable when one considers 
that the effect of a discrete retardation time for a simple Kelvin 
element is significantly felt over only a decade of time. The power 
law provides a continuous distribution of retardation times for the 
retardation spectrum. 

In using the power law for current creep data, the author was 
quite concerned about the extreme sensitivity of the technique to small 
deviations in the strain values. Upon further investigation, the 
power law was found to have a singularity in the vicinity of the 
current data. The acknowledgement of a singularity is quite dis- 
concerting when, from all appearances, the power law is quite well 
behaved. Fig. 3.2 illustrates a typical representation of the variety 
of forms that the power law may assume. In this case, the functions 
are constrained to pass thru (0,0) and (1,1). Fig. 3.3 illustrates 
a different representation of the power law wherein the strains are 






STRAIN 



Fig. 3.3 Power Law Curves Constrained To Go Through Given Strains at Two Non-zero Times. 
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required to pass thru certain values at t-j =1 and t q = 16. These 
times are typical of the time interval span for our experimental creep 
data. A totally different character is now associated with the power 
law. Unlike the representation in Fig. 3.2, the resulting response 
for the case n = 0 is no longer time independent. In fact, this case 
could accurately represent actual creep response over a particular time 
range. 

Equation 3.9 reveals that n = 0 results from z^ ~ ( £ -| + 

For these values of strain, Eqn. 3.10 indicates that e Q is singular, 
as is m calculated by Eqn. 3.11. Because the experimental data falls 
near the singularity, evaluations of the power law parameters are very 
sensitive to small variations in the experimental data. This is 
particularly true when the creep strains are small compared to the 
total strains. 

Figs. 3.4 and 3.5 illustrate the behavior of the power law 
parameters in the neighborhood of the singularity. The figures are 
based on ^ = 1, = 6, t^ = 36, e-| = 1000, and = 1050. These 

are typical values from the current experimental data. At the 
singularity, z^ = 1025, the value of n passes thru zero and the 
values of m and e q diverge without bound. A value of n = .25 was 
typical for the current data. This requires e ~ - 1019.5. The cross- 
hatched regions indicate the errors in computing e Q , n, and m based on 
reading within ± 1/2% accuracy. Because of the disasterous effects 
associated with small errors, precise measurements are imperative to 
the procedure. Nonetheless, the technique is believed to be quite 
useful. As the relative difference between e-j and increases, the 
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Fig. 3.4 Dependence of power law parameters e q and n on for given values of and e^. 
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Fig. 3.5 Dependence of power law parameter m on for given values of e-j and e 
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approach becomes less sensitive to small errors. For specimens 
exhibiting larger amounts of creep, the computed values of the 
parameters were quite consistent. The difficulties arise primarily 
with specimens having a very small transient response. 

The singularity problems do not arise exclusively with the use 
of Eqns. 3.9 - 3.11. Although this approach does contribute to the 
sensitivity, the basic problem is rooted in the nature of the power 
law and how it is employed to fit experimental creep data. The power 
law is not the governing equation for the creep process, but rather 
an empirical technique which has been found to accurately approximate 
the measured creep behavior over a specific time span. The experimental 
data does not suggest any singularities. It is merely that the coef- 
ficients of the power law equation may become singular when attempting 
to model a specific set of data. The power law singularity coincides 
with the transition between a true viscoelastic solid and the region 
labeled as "neither fluid nor solid" in Fig. 3.1. If one assumes that 
a true viscoelastic solid (i.e., infinite creep time results in a 
finite creep strain) exists, it seems obvious that a material could 
also exist with the properties of the transition. Clearly then, the 
singularity is not a result of the technique used to determine the 
power law parameters but is an intrinsic quality of the power law in 
attempting to fit the response characteristics of transition type 
material s . 

No implication is intended that the current epoxy matrix is 
such a transition material. One should realize, however, that the 
power parameters being calculated may be in a steep region of the e Q 
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and m curves in Figs. 3.4 and 3.5. Small errors in reading the experi- 
mental data may result in large errors in the determined values of 
the power law parameters. The power law singularity is significant 
primarily because it explains the sensitivity encountered in deter- 
mining the power law parameters. 

In evaluating the creep data of a specimen subject to several 
stress levels, the values of n, m, and e 0 may be quite inconsistent 
with each other. Realizing, however, that the goal is to express this 
data in the form: 

e(t) - f-j(a) + f-j(cr) t" , (3.12) 

steps can be taken to improve the results. The computed values of m 
may be stabilized by modifying Eqn. 3.11 to remove the dependence on 
the erratic values of e . Knowing that the power law exponent will 
be taken as some constant value n for all stress levels, one may 
write 





(3.13) 


This produces consistent values of m. 

One final comment on the power law is that while it is commonly 

assumed that the value of n is a constant independent of stress, there 

is some evidence that better fits could be obtained by making n a 

function of stress. Fessler and Hyde [29] have found that the creep 

2 

strain of a lead alloy could be expressed up to 14 MN/m by: 

e c = 5.56 x 10“ 6 sinh (a/2.27) t^* 0700 + * 0689a " -00312a 2 ) 
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One set of Griffith's [41] creep data for a 90° specimen indicated a 
linear increase in n with stress. Most of our data, however, was too 
inconsistent to determine functional relationships between n and 
applied stress. There did seem to be a gradual increase in n with 
increasing a, however. 


Principal Orthotropic Properties 

As was noted earlier, the plane stress constitutive properties 
of a linear orthotropic material is completely characterized by four 
independent properties. For the current development, the four 
independent compliances used are: S-j-j, = and Sgg. 

Expressing these as functions of time permits the calculation of any 
viscoelastic response. 

Because the material is nonlinear, however, the response is no 
longer governed by: 


£ i j ( t ’ 6 ) s ijk£ ( t,e ) °k£ (3.14) 

for the creep load 

o ki (t) = a k£ H(t) (3.15) 

where e represents the environmental state such as temperature and/or 
moisture content. For nonlinear elastic materials, the expression for 
strain is often expanded into an odd power series in stress. 

Similarly, for a nonlinear viscoelastic material, one can write 






k£ 


^ijktmnop^’"' 


k£ mn op 


+ 


(3.16) 
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linearization of the compliances about the current stress state, Z . 
This simplifies the equation to 

E ij (t - > * Sijk*' 1 ' ’ E > ^ < 3 - 17 > 

where the compliance tensor has been expressed as a function of the 
stress state. For an orthotropic material under plane stress condi- 
tions, one obtains: 


V,(t) 


S 1 1 ( t , z ) 

S 12 (t,z) 0 


r a. -\ 

°1 
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S 22 (t,z) 0 
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At this point, it should be noted that Yeow [82] and Griffith 
[38] have indicated that the fiber dominated compliances, and 
of the T300/934 material system are neither functions of time nor 
stress, but are linear elastic properties. 

Interestingly, even if the fibers are elastic, as they are 
believed to be, there should be some slight time dependent response for 
S-J-] (and Si ^ ) in a composite material. Sturgeon [73] points out that 
such time dependence of the creep compliance is not due to creep of 
either component, but to relaxation of the matrix in an essentially 
fixed grip configuration. The additional load transferred to the 
fibers results in a small additional strain, which to all outward 
appearances is creep. An important difference between "relaxation 
creep" (as Sturgeon calls it) and creep of the matrix dominated 
compliances is that the former is ultimately limited by the fiber 
response, whereas the latter is dependent on the matrix response and 
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is not limited by the fiber properties. 

While this is an interesting digression, experimental creep 
strains in the fiber direction comprise an imperceptible portion of 
the total strains and may be considered negligible. The elastic 
values of and S^ 2 as obtained by Griffith [38] are used throughout 
the current analysis. Thus it is assumed that only S 22 and S gg will 
be functions of time or stress state. 

The question now becomes what the interaction of stresses will 
be. S 22 , for instance, is clearly a function of a 2 but whether it is 
also a function of a-j or x-j 2 is not apparent. These problems do not 
arise when one considers a uniaxial load state in an isotropic material. 
For combined loadings, or with an orthotropic material, however, inter- 
action effects are an important consideration. Griffith [38] assumed 
there was no interaction effect for the compliances. Thus the 
compliances were expressed as 

s 22 ~ s 22 (t, <? 2 ) 

S 66 = S 66 ^’ T 1 2) 

If one likens creep of nonlinear viscoelastic materials to the 
nonlinear phenomenon of yielding in metals, it would seem appropriate 
to consider interaction effects. In this light, the non-interactive 
expressions of compliance correspond to a maximum stress type yield 
theory. In plasticity theory, the octahedral shear stress, as given 
by. 


t 
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(3.19) 


is commonly assumed to be an accurate simple indicator of the onset of 
yielding. This, of course, is the famous von Mises or Distortional 
Energy Theory of yielding. Thus, in plasticity theory, the octahedral 
shear stress could be considered as the nonlinearizing parameter because 
it controls the onset of the plastic flow. 

The amount of deformation due to the plastic flow is often 
assumed to be governed by the Prandtl -Reuss Flow Rule 


de ij = °i j dx 


(3.20) 


where de.? are the incremental components of the plastic strain tensor 

* \J 

a/ are the components of the deviatoric stress tensor given by 

1 - 1 * 
a. . - -a. • - 6 • . a, , 

ij ij 3 kk 

and dx is an instantaneous proportionality constant. 

Fessler and Hyde [29] note that Johnson, Henderson, and Khan 
have evidence to support an expression of the creep strain as 


£C iO 'f f l<Vt> ^ f 2 (t) ex P(- U/KT ) < 3 ' 21 > 

This expression is consistent with the Prandtl -Reuss Flow Rule and may 
be expressed as 

d ' C ij • "ii ( 3 . 22 ) 

This form is consistent with the incompressibility of creep strains, 
although this assumption is not universally accepted. A recent review 
of creep under combined stresses may be found in Findley, Cho, and 
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Ding [34]. They have noted a paper by Mark and Findley in which it 
was shown that the creep vectors are normal to the Mises ellipse, 
substantiating Eqn. 3.21. 

Lou and Schapery [54] have considered the combined matrix stress 
state induced by uniaxial tests on unidirectional, off-axis specimens 
of glass/epoxy. They have used the matrix octahedral shear stress 
as the nonlinearizing parameter for applying the Schapery integral to 
determine a nonlinear characterization for their material. In their 
work, strain has been expressed as. 
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(3.23) 


Creep in the fiber direction is negligible because the elastic fibers 
are much stiffer than the resin. As a result Eqn. 3.21 is not directly 
applicable to fiber reinforced materials. Furthermore, creep in shear 
[38,54] has been shown to be much larger than that perpendicular to 
the fiber direction. Eqn. 3.23 has been used in the current analysis. 


Matrix Octahedral Shear Stress 

To be a meaningful parameter in predicting the nonlinear aspect 
of the creep of the matrix material, the octahedral shear stress should 
be based on the actual matrix stress. The implication is simply that 
the creep of the matrix should not be dependent on the stresses in 
the fibers. Several approaches have been used to evaluate the matrix 
stresses from the ply stresses. These approaches involve aspects from 
the micro-geometry of the material as well as the individual properties 
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of the matrix and fibers as could be determined from independent 
testing of the two components. 

The simplest approach, known as the mechanics of materials 
approach, models the composite as parallel alternating strips of 
fiber and matrix. This approach results in the rule of mixtures ap- 
proximations for E-j i and 

E n = E f v f + E m (l - v f ) (3.24) 

v 12 = v f V f + v m O " v f) (3.25) 

^22 an< ^ are 9^ ven by 


E f E m 

122 " M 1 - V + E m V f 


G f G m 

’12 - G f (1 - v f ) + G m v f 


(3.26) 

(3.27) 


where 

v^ = fiber volume fraction 

and f and m subscripts denote fiber and matrix properties respectively. 
The matrix stresses based on this model are given by 
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(3.28) 



More exact expressions for averaged values of the matrix stresses 
can be derived, based on more complex modeling of the matrix and 
fiber inclusion. Beckwith [6] has used the Halpin-Tsai equations 
in studying a viscoelastic matrix. Pindera and Herakovitch [57] have 
included residual stress effects and solved for accurate expressions 
for the components of the [B] matrix, based on the use of Hill's [47] 
elasticity approach to reinforced materials. Numerous finite element 
solutions for various fiber shapes and spacings have also been proposed 
(e.g., Foye [36]). 

These techniques all provide improvements over the mechanics of 
materials approach in determining average values of the matrix 
stresses. The major limitation to these techniques is that they are 
based on linear material behavior. For our nonlinear material, the 
appropriateness of averaged values of the ply stresses is questionable. 
The variations of the matrix stresses between the fiber inclusions 
could be quite complex. Indeed, more rigorous expressions for 
averaged values are not necessarily any better for obtaining effective 
values of the matrix stresses in a nonlinear material. 

Besides linear averaging, other averaging schemes can also be 
used. Lou and Schapery [54] have applied a root mean square averaging 
procedure to results from a finite element model. They found this 
approach yielded values of octahedral shear stress which were nearly 
proportional to linearly averaged results for the case of unidirec- 
tional off axis specimens. They also found these values to be in 
fair agreement with the mechanics of materials approach, which was 
subsequently used in their investigations. 
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Independent properties for the matrix and fiber were not 

available. Properties supplied by the manufacturer were inconsistent. 

By assuming a few properties and working backwards from Eqns. 3.24 - 

3.27, the following values were determined and used in the analysis 

E = .487 x 10 6 psi 
m K 

Ef = 34.0 x 10® psi 

v m = - 35 
v f = 63.4% 

E-j ] = 21.7 x 10 6 psi 
= • 33 

While the matrix moduli do decrease with time, they were assumed 
constant in calculating the matrix stresses. 

Once one has obtained the effective matrix stresses, the 
octahedral shear stress in the matrix for a plane stress situation 
can be computed by 

_ 1 v , m rru2 , / rn2 , , itk2 , m,2-,l/2 on N 

T oct “ 3 ^ CT i “ ct 2 ' + ) + \°2 ' + 6 ^ t 12 ' J (3.29) 

Fig. 3.6 illustrates the values of the ply stresses for a uni- 
directional off-axis laminate as a function of fiber angle. Along 
with the typical representations of a-j , and the values of 
a-]™ and T QCt have also been presented. As in Eqn. 3.28, it is noted 
that ^2 (^2 and t *j ^ ^12* 

Adaptations of the Findley Procedure 

The Findley approach to nonlinear viscoelasticity was chosen 
because of ease of application and the available data was sufficient to 
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determine the unknown parameters and seemed to match the data well. 

I 

To evaluate the Findley parameters e Q , o £ , m', and creep 
compliance data is obtained for a given specimen configuration at 
several stress levels. The power law parameters e Q , m, and n are then 
determined for each set of data. 

Findley and his co-workers [31,32,33], and others (e.g.. Boiler 
[7]) have assumed that the value of n is independent of stress level. 
Their published data indicates wide scatter in the value of n. They 
usually take the average value of n to be the constant for the Findley 
equation. Apparently, this approach is based on the inability to 
determine a functional dependence of n on stress because of the data 
scatter. In the present study, the values of n were also widely 
scattered, due in part to the singularity discussed earlier. This was 
a particular problem for short term data where the creep strains were 
small. The values of n were felt to be more stable for specimens 
exhibiting larger creep strains. 

The available data for the unidirectional material, as well as 
the current compliance data for the general laminates, were composed 
of short term (16 or 36 minute) data for several stress levels and 
week long data at some particular stress. Because the long term data 
exhibited larger creep strains, determination of the associated 
Findley parameters was less sensitive to errors than the short term 
data. It was, also felt that by using the value of n from the long 
term data, the model would best predict the overall response of the 
material. Thus the short term data at various stress levels was used 
to evaluate the hyperbolic sine parameters, but the long term data at 
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a single stress level was used to evaluate the exponent n. 

Once a set of e Q and m as functions of the applied stress level 
has been obtained, the hyperbolic sine parameters may be obtained. A 
least squares fit program for a hyperbolic sine function was written 

I 

to provide an accurate and efficient means to obtain e Q , a £ , m', and 
0 . This procedure is described in Appendix A. 

The end result is an equation for strain in terms of the applied 
stress. Knowing the values of the octahedral shear stress at a 
particular load case, one can modify this expression of strain to an 
expression of compliance 

s(t ' T oct> = rf [e 0 s1nh < T oct /0 E > + sinh (WsX 1 < 3 - 30 > 
oct 

In this manner, the uniaxial compliance fitting the experimental data 
at any fiber angle may be determined. 

Determination of Actual Compliance Properties 

Creep compliance data for 90° and 10° specimens at 320°F has 
been obtained by Griffith [41]. This data was used to determine 
expressions for the required and Sgg compliances. Fig. 3.7 
presents the experimental short term compliance for a 90° specimen at 
several creep stress levels. Also included in the figure are the 
results from a week long compliance test at a = 2750 psi . 

The jump in the long term data at t * 100 min was required to 
provide a good fit for the data points. Our original interpretation 
was that there had been a jump in the balance calibration. Conse- 
quently, this strain jump was subtracted from all strain values after 
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t = 100 min. This interpretation has been used to determine the 
power law exponent for the 90° specimen. In light of electrical 
problems encountered in our acquisition of creep strain data, however, 
it was recently decided that poor voltage regulation of the main 
power supply may have caused this discrepancy. 

For most long term creep data from other laminates, the repre- 
sentation of log transient strain vs log time was found to be a fairly 
straight line. This, of course, is in accordance with the power law 
prediction. As seen in Fig. 3.8, however, there is a marked deviation 
from the power law prediction over the time range of 100 to 3000 
minutes. The dashed line in Fig. 3.7 is currently believed to represent 
the best interpretation of the experimental data. The value of n 
predicted from this assumption is 0.22, as opposed to a value of 0.183 
as determined from the strain jump approach. 

Also of concern is that the short time creep rate for the long 
term test is less than half as much as the creep rate for the short 
term data at a comparable stress level. This aspect was not noticed 
until the data were superimposed in the final compilation. The only 
proposed explanation for this disturbing result is that perhaps one of 
the specimens had been postcured and the other had not. An additional 
test or two would be required to correct this anomaly. There is a 
slight difference in the t = 1/2 min comp! i ancestor the long term data 
and the short term data at a comparable load. This is believed to be 
due to small errors in determining the cross-sectional area of the 
specimen, etc., and is not considered significant. 
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As mentioned earlier, the short term data at various stress 
levels was used to evaluate the stress dependence, while the long term 
results were used to calculate the power law exponent. Figure 3.8 
indicates that the data may indeed be plotted along a straight line 
as predicted by the power law. Note that the slope, and hence n, tends 
to increase with increasing stress. Figure 3.9 illustrates a tendency 
to a linear increase in n with increased stress but as mentioned 
before, a constant value was used. Also provided in this figure, are 
plots of e Q and m vs applied stress level. The curves represent the 
least squares fits for hyperbolic sine functions thru the data. The 
fit provided by the Findley approach is considered excellent. It should 
be noted that the transient response is considerably more nonlinear 
than the initial response. 

Figures 3.10, 3.11, and 3.12 represent similar results for the 
10° off-axis data from Griffith [41]. Fig. 3.10 shows that the long 
term compliance is quite similar to that obtained from the short term 
test at a similar stress level. Figure 3.11 indicates that the 
exponent again tends to increase slightly with increasing stress and 
Fig. 3.12 shows that the fit of the Findley approach is again very good. 
The widely scattered values of n are typical of those obtained in 
this study. Published literature indicates a similar degree of scatter 
as may be found in [31,32,33,7]. 

It should be noted that the slopes in Fig. 3.11 vary in a slow, 
regular manner, while the corresponding values of n in Fig. 3.12 are 
very erratic. The explanation is that the values of n plotted in 
Fig. 3.12 are those obtained from Eqn. 3.9. The values of e Q used to 
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determine the transient strains plotted in Fig. 3.11 were the values 
predicted by the hyperbolic sine function which represents the actual 
values of e Q . Had the actual values of e Q , as determined from Eqn. 

3.10 been used, the slopes of the lines in Fig. 3.11 would have been 
the same as the erratic calculated values of n. The use of the pre- 
dicted values of e Q rather than the raw values, results in a moderating 
effect which minimizes the erratic behavior introduced by the 
singularity. 

The Sgg compliance data obtained using the previously discussed 
transformation equations together with the Findley predictions for 
the 90° and 10° data is shown in Fig. 3.13. Because the input data is 
in the hyperbolic sine form, the application of the Findley approach 
to J he calculated Sgg also gives excellent results as shown in Fig. 
3.14. 

Based on the long term compliance data, the corrected value of n 
for the 90° data and the value of n for the 10° data were both equal 
to 0.22. While the agreement for the 10° and 90° exponents is 
excellent, the exponents calculated from Griffith's [38] long term 
30° and 60° compliance data were signi ficantly larger, 0.35 and 0.32 
respectively. If S and S 66 are based on the 90° and 10° results, 
the predicted 30° and 60° response will not be in good agreement with 
the experimental results. The reasons for the discrepancies in 
exponent value have not been investigated. Lou and Schapery [54] have 
indicated that the value of the power law exponent should be 
independent of the fiber angle. It is possible that additional damage 
such as matrix microcracking could have occurred in the 30° and 60° 
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specimens to result in the larger values for n. 

The compliance properties for T300/934 Gr/Ep at 320°F as 
developed in this section may be summarized as, 

= 4.587 x 10" 8 in 2 /! b 
S 12 = -1 .514 x lO -8 in 2 /! b 
S 22 = (l/f)[. 006728 sinh (9/9330) 

+ .00009246 sinh (-9/3030) t’ 1826 ] in 2 /! b 
S 66 = (!/?)[. 009431 sinh (9/8324) 

+ .0002957 sinh (9/2648) f 2162 ] in 2 /! b 
These values have been expressed in terms of the normalized octahedral 
shear stress as given by 

t = 2.410 t . 

oct 

This constant is different from that used in plasticity theory because 
a-j 01 f 0. This normalized octahedral shear stress is equal to the 
applied stress for a 90° specimen, and is used as a convenience. 

Variable Stress State 

The ply stresses are constantly changing with time due to the 
differential creep rates among the plies and possible failure of an 
individual ply. Some type of superposition principle is required to 
account for the strain produced by a variable stress state. The 
Boltzman superposition integral is valid for linear viscoelastic 
materials but not for a nonlinear material. Being an extension of the 
Boltzman integral form to nonlinear materials, the Schapery integral 
expression should directly account for a variable stress history. 

As mentioned earlier, however, this approach was not used because there 
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was insufficient data to evaluate the g o , g-, , g^, and a^ functions. 

The Green-Rivlin approach could be used but, again, is very difficult. 
Furthermore, while these forms do represent a variable stress history, 
their ability to predict strains for varying stresses has not been 
established. 

A variety of simpler approaches have been proposed to account for 
variable stresses applied to nonlinear materials. Most have been 
limited to a uniaxial stress state. Also, they have been used and 
experimentally verified for only one or, in a few cases, several load 
steps. Many of these approaches are graphical in nature so that they 
can account for a general, nonlinear response without requiring a 
characterization of the response. Fessler and Hyde [29] have sum- 
marized several of these basic approaches. 

Fig. 3.15 illustrates the predictions of several graphical pro- 
cedures for a simple step loading given by, 


a(t) - a Q H(t) + (a-j - ct q ) H(t - t-j ) 


(3.31) 


Their application requires knowledge of the independent responses. 


e(t,a 0 ) for a(t) = a Q H(t) 

and 

e(t,a, ) for a(t) = a-j H(t) 

A linear case has been used for illustrative purposes in Fig. 3.15 even 
though the graphical procedures are applicable to any general nonlinear 
response. This permits comparisons with the Boltzman superposition 
integral which is only applicable to linear materials. A linear case 
is schematically represented by letting the a Q response be a constant 
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Fig. 3.15 Comparison of superposition principles for the case of a linear material 
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The time hardening hypothesis requires a vertical (strain) shift 
of e(t,o-| ) to coincide with e(t,a ) at t = t-, . The strain hardening 
hypothesis requires a horizontal (time) shift for the same situation. 
The work hardening hypothesis in general involves a diagonal shift to 
equate the works done in creep. This approach, which is advocated by 
Fessler and Hyde [29], is similar to the strain hardening approach. 

In fact, for a linear material, these two approaches are the same. It 
is interesting to note that none of these graphical techniques simplify 
to the Boltzman integral for a linear viscoelastic material. 

Another approach which has been used for variable stress state 
with nonlinear materials is basically a modification of the Boltzman 
superposition integral . Findley and Khosla [31] have proposed the 
use of a modified superposition principle of the form, 

e(t) = e Q (cSq ) + m(o 0 ) t n 

+ £ 0 ( CT 1 " %) + “ o 0 Ut - t-j ) n 

+ ... 


+ e 0 (°-j - a i_i) + m ( a i " a i_i )(t ■ (3.32) 

for the stress history 

a(t) = a Q H(t) + (a-j - a Q ) H(t - t-j ) + ... 

+ (a. - a i ) H(t - t.) (3.33) 

While good agreement with experimental data was claimed, it 
should be noted that this approach is not consistent. A counterexample 
is found in a simple creep test with 
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e(t) = e Q (a) + m (a) t n (3.34) 

for a(t) = a H(t) 

Now consider a load history of the form 

o(t) = [a/ 2) H(t) + (5/2) H(t - &t) (3.35) 

As At -> 0, the load history approaches that of a simple creep test, 
but the response 

e(t) = e Q ( ct/ 2) + m(a/2) t n 

+ e Q (a/2) + m(a/2)(t - At) n 

= 2[e 0 (S/2) + m(a/2) t n ] (3.36) 

approaches that of a linear material. Obviously this scheme should 
not be used for a nonlinear system. 

In a later paper, Findley and Lai [33] have proposed an alterna- 
tive modified superposition equation of the form 

e(t) = [ £ 0 ( a 0 ) + t n ] 

+ [e 0 <CT., ) + m(o,)(t - tp" - % [a 0 ) - m(o Q )(t - t,)"] 

+ „ . . 

+ Ce 0 (cr i } + Olfaptt - ^ )" - 

- m(a 1 _ 1 )(t - tp"] (3.37) 

While this form does not have the drawback of the previous form, its 
validity for complex stress histories has not been established. None- 
theless, the approach is quite straightforward and easily adaptable to 
a numerical scheme. The method should be used with caution as there 
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are no additional parameters with this approach as exists in the 
Schapery procedure to modify the response for a varying stress level. 
The procedure is a straight superposition of the material response to 
simple creep loading. As such, it requires fewer tests than the 
Schapery integral, but is also less general. 

The author's feeling is that the latter modified superposition 
principle is a simple attempt at predicting the response to a variable 
stress state but its accuracy for complex histories is doubtful. 
Nonetheless, the approach was used in the current investigation. 
Because the variation in ply stresses for creep loading is quite small 
and regular, any valid superposition principle should give reasonable 
results for creep loading. 



Chapter 4 

DELAYED FAILURE MODEL 

Prediction of ply failure within a laminate requires the develop- 
ment of a delayed failure model. Because the state of stress for a 
ply in a laminate is more general than can be modeled by a uniaxial 
specimen, the failure model must account for any arbitrary stress 
state. The approach taken was to generalize an existing orthotropic 
static failure criteria to account for time dependent strengths as 
predicted by existing uniaxial creep strength theories. This combi- 
nation should provide accurate static or very short time failures for 
general loading states, as well as accurate predictions of long term 
failures for a few specific load states. It is then assumed that the 
predictions for long term failures of arbitrary loading will also be 
accurate. 

Most creep rupture criteria for homogeneous isotropic materials 
are based on a linearly decreasing logarithm of the time to rupture 
with increasing stress. This form, as exemplified by the Zhurkov, 
Larson and Miller, and Dorn methods, is given by 

log t r = A - Ba (4.1) 

where t is the time to rupture for a constant creep load of a. A and 
B are material constants for a given temperature [23]. Landel and 
Fedors [50] have noted that in some circles, the form, 

log t r = A - B log a (4.2) 

is viewed more favorably. Because the form of Eqn. 4.2 is a power law, 
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on the creep power law. Furthermore, the relation does not suffer from 
the limitation that a finite rupture time is predicted at zero stress. 
Because of the data scatter and the small range of stresses involved 
with our creep rupture data, however, the preference for one form over 
another becomes academic. Equation 4.1 provides an adequate representa- 
tion for the current data and has been used throughout the present 
analysis. 

Experimentally, the creep stress level is the independent 
variable and the time to rupture at that stress level is the dependent 
variable. For the analysis, however, it is convenient to rearrange 
Eqn. 4.1 to express the creep rupture strength, R, as a function of 
the time to rupture. 

R(t r ) = (A - log t r )/B (4.3) 

The interpretation is that to obtain a failure at time t , one would 
apply a creep stress level as given by 

a = R(t r ) 

Of the several orthotropic static failure theories discussed in 
Chapter 2, the Tsai -Hi 11 criteria was chosen for the current analysis. 
While the Tsai-Wu Tensor Theory is often considered slightly more 
accurate because of its more general form, the requirements of compres- 
sion and biaxial failure data to correctly evaluate the parameters 
prevents wider usage of this technique. While some investigators have 
assumed values for the interaction terms and assumed that tensile and 
compressive strengths are the same, this yields a criteria very 
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similar to the Tsai -Hill approach. Such assumptions do not take 
advantage of the more general nature of the tensor approach and thus 
minimize any advantages over the Tsai-Hill method. 

If the Tsai-Hill criteria is extended to account for time 
dependent creep rupture strengths, the following form results, 

2 2 2 
o' -j 0 -iOq a 0 ii q 

! _ — - J - - — + _ + Lr = ] (4 4) 

[X(t r )] 2 [X(t r >] 2 [Y{t r )] 2 [S(t r )] 2 

Here, the time independent strengths have been replaced by the creep 
rupture strengths which result in failure at t = t . X(t r ) represents 
the creep rupture strength for a uniaxial creep load parallel to the 
fiber direction. For the current material the assumption was made that 
delayed failures do not occur for 0° specimens and that X(t r ) = X. 

Y(t ) represents the functional relation with time of the creep rupture 
strength for a uniaxial creep load perpendicular to the fiber direc- 
tion. S(t ) is a similar functional relation for the shear creep 
rupture strength. Theoretically, S(t ) can be determined from uniaxial 
creep rupture of off-axis specimens and prior knowledge of X and 
Y(t r ). 

Thus, experimental creep rupture data for unidirectional laminates 
is used to determine functional expressions for the creep rupture 
strengths. For creep loading of any arbitrary biaxial stress state 
[ct-j ,a 2 ,Ti 2 ]^ 5 one can compute the associated rupture time by solving 
Eqn. 4.4. 

Griffith [41] has obtained creep rupture data for 90°, 60°, and 
45° off-axis specimens at several temperatures. This data has been 
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replotted with the results for the three orientations at a specific 
temperature on a single graph, and best fit lines have been drawn 
through the data points. The results are given in Figs. 4.1, 4.2, 

4.3, and 4.4. This data was originally presented [38] to show the 
temperature dependence of creep rupture for the various fiber angles. 
The modified rate equation was shown [38] to provide the best 
analytical representation of the data. 

The current interest is in the 320°F data as presented in Fig. 
4.2. The points denoted as "Postcured 60° off-axis" were obtained 
during the present work to determine the effect of postcuring on 
creep rupture. While the magnitude remained about the same, there did 
seem to be a smaller decrease in the creep rupture strength with 
increasing rupture time. Because of the considerable data scatter, 
however, it is not known if this observation is justified. These 
postcured data points were not used for the best fit lines. 

Assuming that the creep rupture strengths may be represented by 
Eqn. 4.3, determination of the slope and intercept for each best fit 
line in Fig. 4.2 allows the determination of the constants A and B. 
Thus one may determine the functional relationship as identified by, 

R 90^V = Y ( t r ) 5 WV and R 45^)- Using the 90 ° creep ru P ture 
strength and one of the off-axis creep rupture strengths, R 0 (t ), one 

should be able to substitute these into Eqn. 4.4 to solve for S(t r ). 

This is accomplished by letting 

o- 1 = cos 2 e R 0 (t r ) 
ct 2 = sin 2 e R 0 (t r ) 
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x-j 2 = sin 0 cos e R 0 (t r ) 

After rearranging the resulting equation, one obtains an expression 
for the shear creep rupture strength, S(t r ), which is not necessarily 
of the form given in Eqn. 4.3. 

This procedure, though straightforward, proved unsatisfactory 
for the available data. Because of considerable scatter in the limited 
amount of creep rupture data available, accurate determination of 
functional expressions for the experimental data are impossible. 

Slight changes in the functions for the 90° and off-axis strengths 
resulted in large variations in the function for the shear strength. 

In fact, if the best fit lines are used for the 60° and 90° data at 
320°F, the predicted shear creep rupture strength increases with in- 
creasing rupture times. Furthermore, the shear strength predicted 
using the 60° data is inconsistent with that obtained using the 45° 
data. Obviously, there are inherent problems in determining an 
independent shear strength relationship by transformations without an 
extensive amount of data. 

To alleviate this situation, the shear creep rupture strength 
was assumed to be of the form 

S(t r ) = a Y(t r ) (4.5) 

Thus, all that must be determined from the data is the value of the 
proportional ity constant, a, as the functional form of S(t r ) has been 
established a priori . There is no rigorous justification to assume 
that the shear strength is proportional to the 90° strength, but such 
appears to be quite reasonable from an intuitive standpoint. 
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Primarily, this procedure reduces the degrees of freedom to a more 
manageable 1 evel . 

To determine the appropriate value of a, the unidirectional 
creep rupture data was again employed. A specific rupture time, t, was 
selected within the range of the available data. Values of the creep 
rupture strengths for this particular rupture time were taken as the 
intercept values of the t = t line and the best fit lines. These 
represent the values of creep stress, for the 90°, 60°, and 45° 
specimens, which would result in rupture at time t. The 60° and 45° 
creep rupture strengths may be normalized with respect to the 90° 
value at that particular rupture time and temperature. 

These normalized creep rupture strengths have been plotted in 
Fig. 4.5 for several times and temperatures. Superimposed upon this 
data are normalized parametric curves representing the Tsai-Hill pre- 
dictions given by Eqn. 2.35 for various values of a. It should be 
noted that these curves will shift slightly depending on the ratio of 
the 90° strength to the 0° strength. For our material, this ratio is 
always very small, and this effect is completely negligible. Based on 
the modified rate equation predictions [38] rather than the best fit 
lines for the creep rupture data, similar results have been presented 
in Fig. 4.6. The octahedral prediction curve in each figure is based 
on failure occurring at a constant value of the octahedral shear 
stress as obtained using equation 3.29. 

These figures are particularly useful in indicating several 
aspects of the failure data. Aside from the parametric curves, the 
points indicate the change in the normalized strengths at a particular 









orientation with respect to temperature. The tendency of the strengths 
at a particular time and temperature to fail along a line of constant 
indicates the appropriateness of the Tsai -Hi 1 1 criteria. A tendency 
for the points at different rupture times to fall along the same curve 
indicates the accuracy of the assumption that the time dependent shear 
strength is a constant proportion of the 90° strength. Presentation 
of information in this form provides a concise yet complete interpreta- 
tion of the data. 

Figs. 4.5 and 4.6 indicate a considerable variation in the 
R 0 (t r )/Rgo(t r ) ratio at the given values of time and temperature. 

Fig. 4.6 illustrates a correlation between increasing a and decreasing 
temperature for the modified rate equation predictions. Furthermore, 
on both figures there appears to be a direct correlation between the 
values of a for and Rgg at a given time and temperature. It is 
of interest to note that the values of a associated with Rgg tend to 
be higher than those for R^g. This represents a deviation from the 
Tsai-Hil 1 equation, which requires the same value of a for all fiber 
angles. Nonetheless, the results for the case of interest, 320°F, 
indicate that the 45° and 60° creep rupture strengths for t = 1 and 
t = 100 minutes are closely clustered around an a = .65 curve. This 
tends to verify that at 320°F the modified Tsai-Hill criteria and 
the assumption that S(t r ) = .65 Y(t r ) are appropriate. Such a formu- 
lation has been used for the failure model. 

It should be noted that the results plotted in Fig. 4.5 (and 
4.6, as well) are still based on drawing a line through only a few 
scattered creep rupture data points. As such, these values are 
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highly subject to error. Verification of the assumed model with any 
degree of confidence would require a more extensive data base and 
is felt to be a priority item for future work in the area. 

The failure properties for T300/934 Gr/Ep at 320°F found from 
the foregoing analysis were as follows: 


X = 195,600 psi 

(4.6a) 

Y ( t r ) = (6800 - 544 log t^) psi 

(4.6b) 

S(t r ) = a Y(t r ) 

(4.6c) 

a = .65 

(4.6d) 


Modification of Failed Ply Stiffness 

The first ply failure does not necessarily result in total 
laminate failure. When failure of one ply has been predicted by the 
failure criteria, the stiffness properties of the ply must be reduced 
to indicate the effect of the failure. Basically, there are two modes 
of failure for parallel fiber reinforced materials. The fibers them- 
selves may break, or the matrix may split along the fiber direction. 

As noted earlier, the Tsai -Hi 11 failure criteria predicts the 
occurrence of failure, but does not predict the failure mode. An 
additional criteria was used to determine the manner in which failure 
occurred. If the failure criteria predicted a failure but the stress 
in the fiber direction did not exceed the 0 ° strength, a matrix failure 
was assumed to have occurred, but no fibers had broken. The matrix 
comp! iances , S 22 and Sgg, were increased by a factor 3 , which was an 
input parameter in the program. Fiber failure was assumed if the 
stress in the fiber direction exceeds the 0° strength. Because so much 



energy is released in fiber failure, the matrix will generally break 
up and total lamina failure results. Thus for a total lamina failure, 
all compliances for that ply were increased by the factor 3. The 
results of the numerical analysis to be presented in Chapter 7 are 
based on g being a large number so that failure effectively eliminates 
the scaled stiffness contributions of the failed ply. Smaller values 
of 3 could be used to allow a failed ply to remain partially effective. 

Cumulative Damage 

Just as a variable stress state causes problems with the non- 
linear constitutive equations, so it complicates the failure analysis. 
While the basic interest at present is to be able to predict creep 
ruptures under a constant load, an important future consideration will 
predict the behavior due to a varying load state. There has been much 
work done for fatigue of composites, but the author was more interested 
in slowly varying loads with only slight variation, as might occur, for 
instance, in a relaxation test. While these predictions for a laminate 
were not intended for the current investigation, they are necessary 
for the ply by ply analysis because of differential creep rates and 
because when a ply fails, the loads are increased in all remaining 
plies. 

A great deal of study has been done in the area of creep rupture 
but nearly all has been based on a constant stress level. The forms of 
the prediction equations do not lend themselves to adaptation for a 
varying stress level. One possibility is the use of some type of 
cumulative damage approach where one takes into account the effect of 
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all the various stress levels acting over their respective interval of 
time. Such approaches have been widely used in fatigue in which it is 
often assumed that failure will occur when 


r n . ^ 


n r-j 


= 1 


where n^ is the number of cycles at stress level ck and n r .j is the 
number of cycles required for a constant amplitude fatigue failure. 
When the parameter m is taken as unity, the familiar Miner's rule for 
fatigue is obtained. Robinson [61] has proposed a similar form for 
damage accumulation in creep, known as Robinson's Life Fraction Rule, 


t • 


= 1 


(4.7) 


where t n - is the time at creep load a. and t r - is the time to creep 
rupture at a ■ . Gerhards [37] has proposed the use of this approach 
for computing residual lifetime in wood. Woo, et al [79] have used 
the approach in predicting cracking in boiler tubes. Davis and 
Coleman [26] have discussed the general conditions under which damage 
may be additive. Kargin and Slonimsky [49] have generalized the 
linear cumulative damage concept to an integral form with temperature 
effect. 


dt 

0 t r ta(t),T(t)] 


1 


(4.8) 


It is of interest to note that Woo, et al [79] have alluded to a 
combined fatigue and creep cumulative damage law of the form 
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Because creep and fatigue damage seem quite different in nature, one 
wonders about the additive nature of the two phenomena. 

Miner's rule for fatigue, as is commonly known, may be grossly 
in error for some applications. The accuracy of a linear cumulative 
damage law for creep is not well established. Preliminary results from 
creep to yield tests on polycarbonate specimens at 75°C indicate that 
this approach may be grossly in error for creep also. There is some 
indication that this is also the case for graphite/epoxy. Primarily, 
the lifetime at a high stress level tends to be increased rather than 
decreased by a previously applied lower load. This is in direct 
contradiction of a cumulative damage theory. This phenomenon is 
believed to be due to some beneficial aging or structural modification 
process. Nonetheless, the linear cumulative damage law for creep has 
been incorporated into the failure model. 

In addition to the cumulative damage law discussed above, a 
number of other approaches have also been proposed. Roberts, Ellis, 
and Bynum [59] have noted several of these approaches, among them: 

Lieberman Strain Fraction Rule 


e H 

Freeman and Voorhees Mixed Rule 






(4.10) 
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and 

Abo El Ata and Finnie Mixed Rule 

e i 

K Z qr — + (1 - K) E — — = 1 (4.12) 

t ri e n 

where e^ is the strain accumulated at ck , e^. is the rupture strain 
at a.. , and K is a weighting constant. Judging from these methods, 
apparently the incorporation of the strain fraction is an important 
addition to the cumulative damage law. While the latter approaches 
are considered more accurate than Robinson's Life Fraction Rule, they 
are also more difficult to use because they require strain data. Use 
of these techniques was not possible as sufficient strain data to 
rupture was not available for the current material. More recent 
developments in damage accumulation are discussed by Bui-Quoc [13] 
and Woodford [80]. 



Chapter 5 

THE NUMERICAL PROCEDURE 

Predictions of time dependent deformations and failures for a 
general laminate require the knowledge of the combined stiffnesses of 
the plies as well as the variation in the ply stresses due to dif- 
ferential creep rates. While the linear elastic properties of a 
general laminate may be obtained quite easily from lamina properties 
by the use of simple lamination theory, there are no simple algebraic 
equations to incorporate viscoelastic behavior. Schapery [66,67] has 
reviewed the use of Alfrey's Correspondence Principle and Laplace 
Transforms to solve boundary value problems for linear viscoelastic 
anisotropic materials. The applicability of this approach to a general 
laminate is based on knowing the time dependent laminate compliance 
tensor. Determination of these functions from compliance properties 
of the individual plies is not trivial. A closed form solution along 
these lines would be very unwieldy to use even if one could be developed. 

DeRuntz and Crossman [27] have developed a numerical procedure 
based on the linear viscoelastic properties of a generalized Maxwell 
element and lamination theory. Crossman and Flaggs [24] have used 
the procedure to analyze time dependent warping of nonsymmetric 
laminates for hygrothermally induced strains. Using the finite 
element method, Foye [36] has modeled the fiber inclusion for linear 
viscoelastic materials. He has averaged the micromechanical response 
to determine the time dependent behavior of a lamina. Using lamination 
theory, this information is incorporated to predict general laminate 
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response. The approach is interesting hut the approximation appears 
to have inconsistent accuracies of the approximations at various stages. 

For the nonlinear model used in the present analysis, the con- 
clusion was made that a numerical scheme incremental in time was 
necessary to predict general laminate response from lamina properties. 
Approaches based on the finite element method and lamination theory 
were considered for this application. 

A finite element approach would have involved discretization of 
each ply into a large number of 3-dimensional elements in order to 
describe the geometry of the test specimen, etc. The advantage of the 
finite element approach would be a complete solution that could predict 
edge effects due to the interlaminar stresses, a 3-dimensional stress 
state distribution throughout the model, and the effects of the end 
constraints. The disadvantages of this powerful approach would be 
the difficulty in developing and using the technique, and the cost of 
solution. To be effective, the grid should be fine enough to reflect 
the variation of the stresses normal to the plane. A refined grid 
would be required all along the free edges and in the vicinity of the 
constrained end. For viscoelastic analysis, a large number of 
storage locations would be required for each element in order to 
reflect the hereditary integral constitutive relations. Use of a grid 
fine enough to show the effect of free edges and constrained ends 
would require such a large amount of storage and computing time, that 
it was undesirable for use at this time. Furthermore, note that the 
stiffness matrix continuously changes with time, which poses a real 
cost problem for the repetitive solution of a large number of 
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simultaneous equations. 

The finite element approach should provide more accurate pre- 
dictions than a lamination theory approach but would represent an 
expensive overkill at the present stage of development. The data 
scatter inherent in creep rupture in general and laminated composite 
materials in particular tends to mask the added accuracy of the FEM. 

Thus while still looking at basic phenomena and trends, the lamina- 
tion theory approach was identified as the most efficient formulation 
to study delayed failures in laminated composites for the present 
effort. 

The Lamination Theory Program 

Classical lamination theory provides a simple means to combine 
the individualized ply stiffnesses to form a unified laminate. As 
mentioned earlier, the fundamental assumption for lamination theory is 
that normals remain straight and normal for any deformation which, in 
effect, causes interlaminar stresses to be neglected. This assumption 
is strictly valid only for interior regions of panels. Without resort- 
ing to other techniques, basic lamination theory results do not lend 
themselves to the analysis around cut-outs or the effect of end 
constraints. In essence, lamination theory provides a technique to 
determine the overall stiffness properties of a laminate based on the 
properties of the constituent laminae. 

A computer program based on lamination theory has been developed 
in the current study. The primary feature of the program is the 
analysis of a general laminate composed of several nonlinear viscoelastic 
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laminae. While the current application has been for constant uniaxial 
creep loads, the procedure was written to handle time-varying biaxial 
loadings as well. Nonlinear viscoelastic laminae properties are input 
internally. The power law has been used in the present work but 
other forms, such as a generalized Kelvin element representation, are 
completely admissible. Compliance interaction effects may be neglected 
or may be assumed to depend on the matrix octahedral shear stress. 

The modified Findley superposition principle (Eqn. 3.37) has been used 
to account for the variable ply stresses caused by differential creep 
rates and in the event of a lamina failure. The modified T sa i -Hill 
failure criteria (Eqn. 4.4) has been used as the ply-wise failure 
theory. The linear damage accumulation rule (Eqn. 4.7) has been used 
to account for the variable ply stresses. All of these procedures are 
subroutines which could easily be modified for other approaches. If 
the initial response is nonlinear in stress, an iterative procedure is 
used to converge on the actual solutions. If a ply "fails", the ply 
properties are modified to reflect the damage state. 

The basic concept behind the incremental lamination theory 
approach is quite straightforward. The solution scheme is based on 
obtaining successive stress, strain, and accumulated damage solutions 
as time is incremented. The strain state is determined at time 
t + At, based on the assumption that the stress state at time t is 
constant for the time step At. New ply stresses are determined at 
t + At based on the current creep strains and the applied mechanical 
load. Accumulated damage in each ply is monitored until ply failure 
is predicted. This procedure is continued until total laminate failure 
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occurs. 

While the actual laminates tested were composed of 16 plies, the 
numerical model need only contain as many plies as the number of dif- 
ferent fiber orientations. This simplification would not be possible 
if one were considering a finite element model or an out-of-plane 
response. 

Numerical Details 

The flow chart in Fig. 5.1 summarizes the numerical procedures 
used in the current approach. The equations for each step can be found 
by reference to Chapter 2. The equations representing the compliance 
and failure models are more involved and have not been stated 
explicitly in the figure. The reader is referred to Chapters 3 and 
4, respectively, for these formulations. The program proceeds as 
follows: Initially, the program reads in the laminae properties, the 

stacking sequence of the laminate, and the control parameters. The 
laminae and laminate creep strains are initialized to zero. The 
instantaneous (hereafter referred to as elastic) laminae compliance 
matrix is determined, based on the current properties and stress state. 
Using lamination theory, the laminae compliance matrices are inverted, 
transformed, and combined to form the laminate stiffness matrix. 

Based on the current applied load, the elastic laminate strains are 
determined. The ply stresses are determined based on the total laminate 
strain and the creep component of the laminae strains. If the ply 
stresses are significantly different than those computed at the 
previous time, a nonlinear iterator is invoked to converge to the 




Fig. 5.1 Flow chart for viscoelastic laminated composite analysis. 
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actual solution. The current stresses are stored in a stress history 
matrix and time is incremented. Assuming ply stresses remained 
constant for the time step, the current ply creep strains are computed 
based on the compliance model discussed in Chapter 3. Determination 
of the total laminate strain requires obtaining the effective laminate 
creep strain. This is determined by first computing the equivalent 
mechanical load necessary to produce the same elastic laminate strain. 
The accumulated damage is evaluated in each ply and compared with the 
failure criteria model developed in Chapter 4. If any ply has 
"failed", the time of failure is calculated and the ply properties are 
modified. This procedure is repeated until all plies have failed. 

Hereditary Integral Evaluation 

Previous efforts in developing numerical procedures for creep 
analysis have primarily been limited to linear systems. For linear 
analysis, the Boltzman hereditary integral 


:(t) = 


D(t - t) 


dg(x) 

dx 


dr 


(5.1) 


has been approximated by a finite series based on discrete time steps 
e(t) = o Q D(t) + (a-, - a Q ) D(t - t-j) + ... 

+ ( a i " a i_-|) °( t “ t-() (5.2) 

This represents an exact solution for a discrete load history given 
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actual solution. The current stresses are stored in a stress history 
matrix and time is incremented. Assuming ply stresses remained 
constant for the time step, the current ply creep strains are computed 
based on the compliance model discussed in Chapter 3. Determination 
of the total laminate strain requires obtaining the effective laminate 
creep strain. This is determined by first computing the equivalent 
mechanical load necessary to produce the same elastic laminate strain. 
The accumulated damage is evaluated in each ply and compared with the 
failure criteria model developed in Chapter 4. If any ply has 
"failed", the time of failure is calculated and the ply properties are 
modified. This procedure is repeated until all plies have failed. 

Hereditary Integral Evaluation 

Previous efforts in developing numerical procedures for creep 
analysis have primarily been limited to linear systems. For linear 
analysis, the Boltzman hereditary integral 


e(t) 





(5.1) 


has been approximated by a finite series based on discrete time steps 

e(t) = a Q D(t) + (a-j - o Q ) D(t - t ] ) + . . . 

+ (a. - a i _ 1 ) D(t - t.) (5.2) 

This represents an exact solution for a discrete load history given 
by 
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o(t) - cr Q H(t) + (a^ - a Q ) H(t - t, ) + ... 

+ (a i - a i _ ] ) H(t - t.) (5.3) 

As indicated in Chapter 3, the procedure used in the present approach 
is based on the modified Findley superposition principle of the form 

e(t) = a Q D(t,G Q ) 

+ G 1 D(t - tpa-j) - c Q D(t - t-j ,a Q ) 

+ ... 

+ ^ i D(t - t • »a . ) - a._ 1 D(t - t^^^) (5.4) 

Equation (5.4), while not expressed as a convolution integral, is con- 
sidered to be a discretization of a modified hereditary integral. 

Thus, subsequent solution techniques are referred to as a numerical 
evaluation of a convolution or hereditary integral. 

Convolution integral evaluations require large amounts of 
storage locations to keep track of the time history effects. For 
linear materials, Zienkiewicz [86,87] has proposed the use of a 
generalized Kelvin model for the viscoelastic behavior. One may 
express the creep strain increments in terms of the current stress and 
the individual creep strains for each Kelvin element 



[ 2 L.’ 
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i-i n i c j 


where Ae c is the increment of creep strain 

■f" h 

E. - spring modulus of i Kelvin element 

j. i- 

n.j - dashpot coefficient of i L Kelvin element 
a - applied stress 


(5.5) 
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e^, 1 - creep strain in i^ Kelvin element 
and N - number of Kelvin elements. 

This procedure is more computationally efficient than a direct 
integration of the convolution integral. Taylor, et al [76] have 
extended this approach to account for changes in the material properties 
due to temperature, etc. One problem with the technique is that of 
fitting a material with a generalized Kelvin model. To express an 
increasing compliance over many decades of time requires a large 
number of Kelvin elements to provide a good fit. Crossman and Flaggs 
[24] have shown a very wavy representation for the compliance of an 
epoxy matrix using 10 Kelvin elements over a 20 decade time span. If 
10 to 20 Kelvin elements for several compliance terms of each ply are 
necessary, a large amount of storage is required as well. This 
approach was originally used in the current investigation but the book- 
keeping became difficult. Keeping track of all the individual Kelvin 
element properties and creep strains was cumbersome. The main limita- 
tion of this approach, however, was that it is only valid for a 
linear material. Therefore, this procedure was eventually abandoned 
because it could not model the nonlinear behavior of the epoxy matrix. 

The approach actually used, then, was a direct numerical integra- 
tion of the modified convolution integral (if it could actually be 
expressed) by the finite sum given in Eqn. 5.4. By using logarithmic 
increments in time, the number of stress history storage locations 
required was quite reasonable. 



105 


Iterative Scheme for Nonlinear instantaneous Response 

If e Q is not linear in a, the material system exhibits non- 
linear instantaneous response. As such, an iterative scheme must be 
used to determine the ply stresses based on the applied load. At 
first application of the load, ply stresses are calculated based on an 
arbitrary value of x. Based on the ply stresses predicted from this 
first computation, revised values of x are calculated and the procedure 
is repeated until convergence is achieved. Convergence is checked by 
a square norm 
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( 5 . 6 ) 


Because the material is only slightly nonlinear in the 
instantaneous response, convergence is achieved quickly and easily by 
this scheme. A more sophisticated scheme is required for highly non- 
linear materials. This iterative convergence scheme is also required 
if the applied load changes, or if a redistribution of ply stresses 
occurs because of a ply failure. 


Log Increments in Time 

Because of the nature of creep response, the creep rate decreases 
significantly as time increases. Logarithmic rather than real time 
increments are used. Based on the rate of activity, this applies a 
more consistent amount of attention to the time regions. Thus, by 
running ten iterations per decade of time, for example, a more 
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consistent base is used for the time increments. This permits close 
scrutiny near the initial loading and minimizes computational efforts 
when analyzing the long time regions where the creep rate is extremely 
slow. 

When a ply fails, the creep rates immediately following may be 
quite high as the stresses in each remaining ply adjust to carry the 
applied load. Thus, when a ply fails, the time step is reduced to 
carefully follow the action after failure of a ply. The logarithmic 
incrementation scheme is again employed, but with a constant time added 
to it equal to the time of ply failure. 

Fiber Rotation Due to Large Deformations 

Some specimen configurations may exhibit considerable deformations 
and corresponding fiber rotations. As the specimen elongates in the 
axial direction and contracts in the transverse direction due to the 
Poisson effect, fiber angles tend to become smaller, shifting more of 
the load to the fibers, effectively making the specimen stiffer. To 
account for this effect, one can calculate the amount of fiber rotation 
and use the new fiber direction in subsequent calculations. Assuming 
a general deformation of an element containing fibers at an arbitrary 
angle e to the x axis, it can be shown that the new fiber angle 6' is 
given by: 


e 1 = tan 


(tan 9 ) (1 + e y ) 

1 + E x + Y xy tan 8 


( 5 . 7 ) 


This procedure is applied to all plies at each time step. 
Assuming that the old fiber angle can be replaced by the new angle 
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represents a first order approximation of the actual effect. It 
should be noted that 0 in Eqn. 5.7 must always be taken as the ply 
angle of the original layup. Use of the updated values would produce 
incorrect results. 

Residual Thermal Stresses 

Because of the considerable difference in the thermal expansion 
coefficients in the fiber and matrix directions, significant residual 
thermal stresses in general laminates may be induced by the cool-down 
process of the cure cycle. In fact, linear elastic analysis of this 
aspect may predict ply failure during cool-down and prior to the 
application of mechanical loads. Consideration of residual thermal 
stresses is often felt to be important in laminate strength predictions 
but have not been considered in the current numerical scheme. As dis- 
cussed in the development, thermal stresses may be included by adding 
the independent thermal strains to the creep strains in Eqn. 2.9. All 
of the current work was done at elevated temperatures, and most were 
done at 50°F (28°C) below the glass transition temperature. At these 
high temperatures, the residual stresses should be quite small. 
Furthermore, Weitsman [78] has shown that a viscoelastic analysis yields 
residual stresses about 20% smaller than those predicted by linear 
elastic analysis. The soak at elevated test temperatures prior to 
testing should lead to even further relaxation of the residual stresses. 
Thus the thermal stresses have been assumed to be negligible for the 
current work, but they may need to be considered for work at lower 
temperatures . 
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Numerical Stability 

As with many iterative numerical schemes, instability problems 
were encountered with the program. A minor problem arose from 
numerical oscillations produced by the sudden application of the 
initial mechanical load. Typical disturbances of the ply stresses were 
on the order of 10% and were damped out within about 8 iterations. 

This problem was improved by using a technique somewhat analogous to 
the central difference method. The octahedral shear stress in each 
ply was moderated by 

T oct = C*Oold + (x oct ) new ]/(1 + ^ (5 ’ 8) 

Letting 4) = 1 results in a straight average of the previous and the 
current octahedral shear stress values. The "old" value conies from the 
previous iteration and the "new" value represents the value based on 
the current ply stresses. The result from Eqn. 5.9 is an averaged 
value which improves stability considerably. 

The major instability problem occurred only when analyzing the 
two fiber orientation laminates. Theoretically , the procedure should 
be stable as long as the time steps are sufficiently small. The 
computed values of the octahedral shear stress in the two ply laminates 
are larger than those encountered in the three or more ply laminates. 
The nonlinear dependence is based on a hyperbolic sine function in- 
volving -r oct . As T Qct becomes larger, relatively small changes in 
T oct resu ^ in huge variations in the hyperbolic function. This pro- 
vides the driving force for the instability. As seen in Figs. 3.7 
and 3.10, the compliance turns up sharply at longer times when plotted 
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against log time. Thus with the use of logarithmic time increments, 
the compliance curve becomes harder to follow accurately at long 
times. This results in the initial perturbation of the solution. 

Increasing if) in Eqn. 5.9 to values on the order of 100 produces 
stable results but leads to a very sluggish response. However, when 
proper time increments are chosen, stable results are obtained over 
long time spans for most laminates. 



Chapter 6 

EXPERIMENTAL PROCEDURES 

Several phases of experimental work were conducted so that the 
phenomena of creep and creep rupture in laminated composite materials 
would be better understood. The main emphasis was on obtaining creep 
rupture data for several different layups which heightened our under- 
standing of creep rupture of general laminates and helped identify the 
significance of time dependent behavior for practical applications of 
laminated composites. Data was also furnished for comparison with the 
program predictions. Short and long term creep compliance data was 
taken at several stress levels for several specimens. This provided 
information on the nonlinear tendencies of general laminates, as well 
as a check on the program results. 

The Material and Specimens 

The material used for the current experimental work was Union 
Carbide T-300 graphite fibers preimpregnated with Fiberite 934 epoxy 
resin. The 16 ply panels were manufactured by Lockheed Corporation 
and had a nominal thickness of .086". Most of the material was pro- 
cured in 1979 and will be referred to as "new" material. One panel 
was procured in 1978 and will be referred to as the "old" material. 

This is consistent with the notation used by Griffith [38], Unless 
otherwise noted, the material is assumed to be from the new batch. 

The specimens used for creep compliance and creep rupture testing 
were sawed from the panels with a diamond abrasive disk. Nominal 
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specimen dimensions were 10" x 1 / 2 " which resulted in a gage length of 
6-1/2" outside the grips. This produced an aspect ratio of 13 
which is normally considered adequate for testing orthotropic materials. 
Alignment holes (1/8") were drilled in each end of the specimens. 

A variety of specimen layups was desired to provide a wide base 
of experimental data. Because of the expense and time lag in obtaining 
additional graphite/epoxy panels, the specimens were cut at several 
off-axis angles from the panels on hand. While this permits many 
specimen configurations, the off-axis specimens are in general un- 
balanced about the test direction. The problems associated with this 
aspect are discussed in the following chapter. 


Specimen Configurations 

All specimens used in the current work were cut at various angles 
from one of the four available parent panels: 


Panel 

#1 

[O/30/-3O/O] 2s 

"old" 

material 

Panel 

#2 

Was 

"new" 

material 

Panel 

#3 

[0/45/-45/0] 2s 

"new" 

material 

Panel 

#4 

[0/90] 4s 

"new" 

material 


While panel #4 has been referred to throughout this paper as 
[0/90]^ s , it was recently found that the actual layup is 
[0/90/0/90/90/0/90/0] s . The different specimen configurations were 
each assigned a designating letter as follows: 

A - [90/60/-60/90] 2s 

These specimens were cut at 90° from panel #1. These speci- 
mens are differentiated from the C specimens noted below 
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because the scrim cloth was not removed prior to the post- 
cure. The post-cure deteriorated the scrim cloth on each 
side of the specimen to a point where it would not peel 
off. It could only be removed by tediously flaking off 
small pieces. It was decided to leave the deteriorated 
scrim cloth on the specimens rather than risk damaging the 
specimen surface. The results from these specimens were 
considered equivalent to the C specimens. 

B - [60]g s cut at 60° from panel #2. 

C - [90/60/-60/90] 2s cut at 90° from panel #1 . 

D - [75/45/-75/75] 2s cut at 75° from panel #1 . 

E - [10/55/-35/10] 2s cut at 10° from panel #3. 

F - [20/65/-25/20] 2s cut at 20° from panel #3. 

G - [90/45/-45/90] 2s cut at 90° from panel #3. 

H - [75/30/-60/75 ] 2s cut at 75° from panel #3. 

I - [60/1 5/-75/60] 2s cut at 60° from panel #3. 

J - [15/-75]^ s cut at 15° from panel #4. 

K - [30/-60]^ s cut at 30° from panel #4. 

Panel #3 is similar to the layups used in current applications 
of graphite/epoxy to military aircraft wing skins, etc. Specimens E 
and F represent the situation of the principal stress being slightly 
off axis to the principal axes of the material. Such loading is 
particularly pertinent to the stress distributions around cutouts and 
pylons, and near leading and trailing edge fittings. 
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Post-Cure 

In studying the "old" material, Yeow [82] found no dependence 
of the material response on the thermal conditioning, or aging, of 
the specimens by holding them at an elevated temperature for given 
periods of time. Griffith [38] however, found that for his specimens 
cut from the "new" material, there was significant dependence of the 
material compliance on the amount of time the specimen was held in the 
test oven before testing. He concluded that there was a significant 
post-curing effect occurring in the "new" material. It was also dis- 
covered that the "old" batch of material had been subjected to a post- 
cure process, but that this extra step had been eliminated from the 
manufacturing process for the "new" batch of material. 

In an effort to reduce this aging effect of the specimens while 
under creep loading in the oven and to insure that previous condition- 
ing of all specimens was similar, all specimens used in the current 
work were subjected to a thermal conditioning cycle shortly after 
being cut from their parent panel. The post-cure cycle was chosen to 
give a gradual heat-up period from room temperature to 380°F at a 
rate of 50°F/hour, followed by a 36 hour hold at 380°F, and then a 
very slow cool-down to room temperature at a rate of 5°F/hr. This 
very gradual cool-down helped minimize the lock-in of excess free 
volume and reduced the residual stresses. 

A cam was cut for the Blue H oven to provide the required 
temperature for the cycle. The specimens, with end holes already 
drilled, were hung vertically on wire hangers in the oven. The 
hangers were made in such a way to space the specimens apart to prevent 
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them from colliding because of the forced convection currents of the 
oven. Many of the specimens became bowed as a result of the post- 
cure cycle. This warping was limited to specimens cut with matrix 
dominated test directions. Specimens of a given configuration all 
warped in the same direction and about the same amount. Fiberite 
indicated that this warping may result when the plies are not nested. 
Apparently, a ply is not symmetric about its midplane, but has a pre- 
ferred side and if all the plies of a panel are stacked with the same 
side towards the tool, the panel will also not be symmetric about its 
midplane and could result in warping in the post-cure process. While 
the warping was not considered particularly detrimental for the 
specimens, it is thought that this effect can be eliminated by properly 
nesting the plies in the layup process. Nesting implies alternately 
placing the backing side of the plies toward and away from the layup 
tool . 

Moisture Content 

Moisture content is known to have a significant effect on the 
strength and stiffness of the epoxy resin. Absorption of moisture is 
believed to have an effect similar to lowering the glass transition 
temperature [45]. Thus, in light of the Time Temperature Super- 
position Theory, moisture content is seen to be an accelerating 
parameter for the creep process. Obviously, therefore, moisture 
content could be a factor in creep studies but was not studied directly 
in the current investigation. To minimize the effect of moisture, 
efforts were taken to maintain consistent moisture levels in the 
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specimens . The specimens were placed in 3 desiccator upon removal 
from the post-cure oven and stored there until tested. The moisture 
content of the air within the desiccator was maintained at 13-18% 
relative humidity. 

The moisture content of several specimens was monitored over a 
6 month period. The specimens lost about 1% of their total weight 
during the post-cure cycle, presumably moisture. There seemed to be 
a very slight decrease in weight (.05%) over the first few days in the 
desiccator. The only explanation for this phenomena is that the long 
cool-down period for the post-cure cycle permitted humid room air to 
enter the oven and provide moisture for the specimens to absorb. The 
desiccator air was drier and removed some of the moisture. The weight 
of the specimens then remained essentially constant throughout the 
6 month period. Thus, moisture content of all specimens was felt to 
be quite consistent. 


Equipment 

Three lever arm creep machines were used for creep loading of 
the specimens: a 6000# Budd machine, a homemade frame similar to the 

Budd, and a 20,000#, Series 2330, ATS machine with automatic releveler. 
Each machine was equipped with an ATS Series 2912 oven and Series 230 
temperature controller. These control units were very temperamental 
but when working, maintained oven temperatures within ± 2°F (± 1°C) 
of the desired temperature. Temperatures were determined with a 
Doric Model 41 2A Trendicator. Hewlett Packard 7100B Strip Chart 
Recorders were used to record the strain data for the creep compliance 
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testing. A 2100 Vishay System was used to condition the strain gage 
output. Shunt resistors were used for calibration of the strain 
gages. 

Fluctuations in line voltage plagued the collection of creep 
compliance data. The Vishay conditioning unit could not provide 
adequate compensation for the erratic line voltage. The 110V line in 
the lab surged in the evenings and on weekends as if from reduced 
power consumption during off hours. To minimize this problem, a SOLA 
Constant Voltage Transformer was connected to the 220V line. The off 
hour power surges were eliminated with this arrangement, but large 
voltage spikes at one minute intervals would sometimes occur on this 
circuit. These spikes were very annoying, but did not disturb the 
data as much as the power surges. This phenomena was apparently due 
to the operation of some piece of equipment which was operational 
primarily on Friday through Monday. Several other experimenters in 
the building indicated that they have experienced similar phenomena. 
While the source of the problems could not be located, it is recom- 
mended that the situation be remedied before further creep testing 
is undertaken. 

Creep Rupture Data 

Before testing, the specimens were measured to determine the 
cross-sectional area used to compute the stress level. Failure was 
assumed to occur at the location of minimum cross-sectional area within 
the specimen length. The minimum cross-sectional area may not cor- 
respond with either the minimum width or thickness when both vary 
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over the length of the specimen. Because this evaluation of the 
minimum cross-sectional area would be very tedious, it was assumed 
that one could scan the length of the specimen for the minimum thick- 
ness (width), measure the width (thickness) at that location, and use 
these values to compute the minimum cross-sectional area. 

The measured specimen was mounted in clamping grips and placed 
in the test frame. At least one half hour was allowed before the test 
was begun so the specimen could reach thermal equilibrium at the test 
temperature. Full application of the load was spread over several 
seconds so that the dynamic effects were minimized. The ATS machine 
at the 20:1 lever arm ratio tended to jerk the specimen at loading 
and during lever arm releveling. This caused oscillation of the load 
train system. This was of some concern, but was found to be 
negligible if care was used in loading the specimen. See Appendix B. 

Previous efforts indicated that end tabs bonded on unidirectional 
and [90/60/-60/90] 2 $ specimens precipitated creep rupture at the end 
tab. Based on these findings, end tabs were not used on most creep 
rupture specimens. However, without end tabs, the two strongest con- 
figurations (specimens E and F), however, failed through the hole 
within the grips. To prevent such failures, cross ply glass/epoxy end 
tabs were bonded to these specimens with Epoxylite 5403 (M-Bond GA-61 ) . 
This adhesive requires a high temperature cure cycle. To prevent cool 
down anomalies, the post cure cycle was repeated but with a 3 hr hold 
at 380°F. The tabs were successful in preventing grip failures in 


most cases. 
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Specimen Measurement Difficulties 

Probably the most difficult part of creep rupture experiments is 
the measurement of the cross-sectional dimensions of the specimens. 

The problem is to measure the specimen at the precise location before 
loading where fracture will occur after loading. Once the failure has 
occurred, the specimen cannot be reliably measured because the delamina- 
tions and splintering may extend for a considerable distance beyond 
the actual separation. Thickness measurements within these debonded 
regions would be meaningless. If one goes far enough away from these 
zones to obtain an accurate thickness measurement, the thickness at 
the new location may be considerably different than that at the 
fracture point. 

Part of the problem lies in the existence of thin spots 
occurring in many of our specimens. These regions were, in some 
cases, 1/2" long over which the thickness tapered so as to be .003"- 
.005" thinner than the nominal thickness. Such thin spots were 
believed to correspond to splice locations where the edges of two 
plies were not quite brought together. During the cure process, resin 
flowed in to fill the void. If all the splices on several ply levels 
occurred at the same locations, as in mortar joints in a brick wall, 
the resin flow could result in a significant thin spot. For some 
specimens, the thin spots were noticeable enough to be located merely 
by passing the fingers over the specimen. When splice lines were far 
enough apart, it was sometimes possible to cut the specimens in such 
a manner that the thin spots did not occur within the test length of 
the specimen. If, however, one is trying to get as many specimens as 
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possible out of a panel at some angle to the splice direction, that 
is not always an option. The best solution to these problems in the 
future would be to persuade the manufacturers to stagger their 
splices, in spite of possible claims of increased waste. 

Often in composite work, the specimen strength is thought to 
be somewhat insensitive to the thickness. In fact, in some circles, 
the thickness is often assumed to be some nominal value for all 
specimens. The concept is based on the assumption that a thick or 
thin spot indicates only the amount of matrix present while the fiber 
content remains fairly constant. This is probably a good assumption 
for specimens with fiber dominated strengths, but not for matrix 
dominated specimens. 

Specimens do not always fail at the thin spots, because the 
actual load situation is considerably different than the idealized 
creep test of a uniform homogeneous material under a uniform uniaxial 
stress state. As stated previously, for a specimen of varying width 
and thickness, the smallest corss-sectional area does not necessarily 
occur at the location of either the minimum width or thickness. 
Further, the specimen contains flaws and defects which may render the 
specimen weakest at still some other location, than the one for 
minimum area. Also, particularly in the unbalanced specimens, the 
stress state is neither uniform nor uniaxial. Furthermore, especially 
for specimens with small fiber angles, it is possible for a failure 
initiating at one location to result in a final separation some 
distance away. These reasons indicate that the strength of laminated 
specimens may not be as sensitive to cross-sectional area as a 
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homogeneous, isotropic material might be. This could result in 
greater data scatter for composite specimens than is observed for 
metals or polymers. 

Crooked Fibers 

An interesting anomaly was noticed in many of the specimens cut 
from the "old" [0/307-30/0^ panel. Many specimens contained fibers 
with a large degree of curvature in certain regions of certain plies. 
These plies exhibited deviations from the desired fiber angle by as 
much as 30°. These sharp bends occurred over a very local region, but 
could indeed have had a significant effect on the strength of a 
specimen. Typical photographs of some of these specimens are shown 
in Plates 6.1 and 6.2. An examination of all the specimens from this 
panel which were creep ruptured showed that those specimens exhibiting 
crooked fibers at the fracture zone had often failed at shorter than 
predicted time to rupture. While the correlation was not perfect, it 
was very significant. Obviously, the stress state induced near these 
knees could be considerably worse than the nominal stress. Also, if 
the fiber angle increased, the laminate would lose part of its load- 
carrying capacity. It should be noted that some of the specimens 
which broke prematurely did not externally exhibit significantly 
crooked fibers. However, only a few plies are actually visible at the 
failure zone. Possibly, therefore, some of the internal plies were 
crooked and precipitated the failure mechanism. 

There was little that could be done to account for the effects of 
crooked fibers. What portion of the manufacturing process of the 




Plate 6.1 Crooked Fibers. 



Plate 6.2 Crooked Fibers. 
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[0/ 30/ -30/0^25 panels which caused this phenomena is not understood. 
Crooked fibers were not noticed in other panels. Our recommendation 
is for more careful processing of all composite panels. If com- 
mercial producers cannot provide better products, it may be desirable 
to obtain the equipment to produce panels in the department under 
more careful supervision. One would still wonder about the quality of 
products used in practical applications, however. 

Baseline Data for Creep Rupture 

Baseline strength data is needed prior to taking actual creep 
data in order to minimize the trial and error procedure necessary to 
determine the stress range for which delayed failures will occur in a 
reasonable amount of time. Baseline static strengths were determined 
from constant crosshead tests on an Instron Test Machine. Fig. 6.1 
illustrates the baseline data for a [90/60/ -60/90l 2s laminate taken 
at several temperatures. Note that the strength declined as the 
temperature increased. 

Creep Yield of Polycarbonate 

Because Gr/Ep is expensive, difficult to obtain, and hard to 
machine, finding other materials which can be used to investigate 
basic rate processes in general is desirable. Polycarbonate is in- 
expensive, readily obtained, and easily machined and had been pre- 
viously studied in our laboratory by Brinson [8]. While brittle, 
fiber reinforced epoxy seems quite different than ductile, isotropic 
polycarbonate, there are several similarities in their time dependent 
response. Therefore, polycarbonate dogbone specimens have been tested 
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Fig. 6.1 Ultimate stresses for laminate A at constant crosshead rate for various 
temperatures . 
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to obtain creep yield failure data. Yield was defined as large scale 
yielding or Luder's band formation as measured by an extensometer and 
strip chart. The results from this work along with the repercussions 
for Gr/Ep are discussed in the next chapter. 

Strain Measurements 

Griffith [38] investigated several methods of preparing Gr/Ep 
specimens for mounting strain gages to be used for determining creep 
compliance and the procedure he described was used to mount the strain 
gages for our compliance testing. Two 350ft gages were mounted on 
opposite sides of each specimen and wired in series. This produced an 
effective 700ft configuration which minimized gage heating and elimi- 
nated specimen bending effects. Because only one specimen from each 
laminate was strain gaged, a specimen from the same laminate but not 
the same fiber orientation was used for thermal compensation. The 
errors introduced by this aspect were completely negligible for the 
± 2°F temperature variation maintained during the compliance testing. 



Chapter 7 

RESULTS AND COMPARISONS 

This chapter contains the experimental and numerical results of 
the current investigation. The creep rupture data is reported first 
to permit an independent assessment of the delayed failure results. 

Creep compliance data--both experimental and as predicted by the 
lamination program--are presented and compared. Predicted creep 
rupture strengths are then presented and compared with the experimental 
data. Photographs of typical delayed failure zones are shown along 
with photomicrographs of edgeviews of several specimens. The poly- 
carbonate results are indicated and similarities in the rate processes 
and physical aging effects of polycarbonate and Gr/Ep are noted. 

Comments about the accuracy of predictions are then given. Finally, 
a discussion of the grip constraint stresses in the unbalanced laminates 
is presented. 

Creep Rupture Data 

Obtaining creep rupture data for several laminates was a primary 
thrust of the experimental phase of this work. One hundred-ninety 
specimens from nine different laminates were prepared and tested to 
obtain creep rupture data. Much of the data shows considerable 
scatter, as is typical of creep rupture data in general. Nonetheless, 
delayed failures were produced in all of the nine different laminates 
tested, substantiating the claim that creep and delayed failures are 
of real engineering concern. 
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The creep rupture data for each laminate is presented in Figs. 

7.1 - 7.12. The figure legends indicate the type of failure that 
occurred. "Good failure" indicates that the separation was located 
well within the test region. "Failed near grip" signifies that the 
fracture occurred within the test region of the specimen, but was 
very near the grip. While most of these data points are believed to 
be valid, it should be noted that the stress state near the grips may 
be more complex. This is particularly true for the unbalanced laminates 
in which the stress state induced by the grip constraints may be 
considerably different than the desired uniaxial stress state. 

"Failed in grip" denotes a failure within the grip. These data points 
are not viewed as being representative creep ruptures. The "+" symbol 
denotes an A specimen tested with the scrim cloth intact. These 
points are comparable with the C specimen results and they have been 
plotted together. 

All pertinent data points have been depicted in the figures. A 
least squares fit (LSF) straight line has been drawn through the data 
points which are believed to be valid. Those points not considered as 
representative creep ruptures have not been used in the least squares 
fit and have been darkened in on the figures. Such points include: 

1) failure within a very short time of loading (20 sec.) 

2) failure within the grip 

3) specimens with obvious crooked fibers (laminates A and C 
only) 

4) specimens which did not fail 

5) specimens with loading anomalies. 
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Fig. 7.1 Creep ruptures of A and C specimens ([90/60/-60/90] ? ) at 290°F (143°C) 
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Fig. 7.2 Creep' ruptures of A and C specimens ([90/60/-60/90] 2s ) at 320°F (160°C). 


A 


A FAILED NEAR GRIP 
O GOOD FAILURE 

FILLED IN POINT NOT 



CREEP RUPTURE A,C SPECIMENS AT 320°F 

(I60°C) 

[90/60/-60/ 90 ] 


2s 


62 


58 


54 


50 


STRESS (MPa) 



STRESS ( ksi) 



62 


58 


54 


50 


LOG RUPTURE TIME (min) 


Fig. 7.3 Creep ruptures of A and C specimens ( [90/60/ -60/90]2 S ) at 350°F (177°C). 
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Fig. 7.6 Creep ruptures of E specimens ([10/55/-35/10] 2s ) at 320°F (160°C). 
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Fig. 7.7 Creep ruptures of F specimens ( [20/65/-25/20] 2s ) at 320°F (160°C). 
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Fig. 7.8 Creep ruptures of G specimens ( [90/45/45/90^ ) at 320°F (160°C). 
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Fig. 7.9 Creep ruptures of H specimens ([75/30/-60/75] 2s ) at 320°F (160°C). 
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Fig. 7.10 Creep ruptures of I specimens ([60/1 5/-75/60] 2$ ) at 320°F (160°C). 
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Fig. 7.11 Creep ruptures of J specimens ([-15/75]^) at 320°F (160°C). 
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. 7.12 Creep ruptures of K specimens ([-30/60]*,,) at 320°F (160°C) 
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The cross-sectional areas for specimen configurations A, B, and 
C were based on the dimensions at the narrowest width. It was later 
realized that this was probably the wrong measurement method for these 
matrix dominated laminates. Because of the large fiber angles, the 
damage zone was quite short and the dimensions near the fracture could 
be obtained. Creep rupture data based on these failure dimensions was 
plotted in addition to data based on pre-test dimensions. While 
individual data points were slightly shifted, the best fit line and 
degree of data scatter were found to be similar for the two measurement 
techniques. Only the data based on the original dimensions is 
presented herein. 

The cross-sectional areas for specimen configurations D, E, F, G, 
H, I, J, and K were based on the dimensions at the smallest thickness 
within the test length. A dial indicator with .0001" sensitivity was 
used to facilitate scanning the lengths of these specimens for the 
thinnest location. Watching the dial while passing each specimen 
under the indicator stylus provided an efficient means to locate the 
minimum thickness. 

Creep rupture data based on a nominal thickness was also plotted 
for each specimen. This was, in general, found to be less satis- 
factory than the creep rupture data based on either the initial dimen- 
sions or the failure dimensions. For a general laminate, the effective 
thickness is probably somewhere between the actual thickness and some 
nominal value, as was discussed in Chapter 6. 

A and C specimens were tested at three different test tempera- 
tures. Figs. 7.1 - 7.3 illustrate creep rupture of the A and C 
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specimens at 290°, 320°, and 350°F, respectively. In Fig. 7.4, the 
valid creep rupture data points from the three previous figures have 
been superimposed along with their best fit lines. It is observed 
that the lines are parallel and nearly equally spaced. The modified 
rate equation (Eqn. 2.41) was applied to these lines and very good 
agreement was obtained as expected when the three parameter equation 
is used. 

All subsequent creep rupture data was taken at 320 C F. This 
temperature provides a good compromise for staying well below the 
glass transition temperature yet above the temperatures at which the 
creep compliance and creep rupture curves flatten out. Fig. 7.5 
illustrates the creep rupture of the D specimens. Data for the E 
laminate is found in Fig. 7.6. As noted earlier, the latter specimen 
represents a principal stress 10° off axis from the primary fiber 
orientation of a practical laminate. Delayed failures were obtained 
and indicate a noticeable decrease in strength with time for creep 
loading. The results from the F specimens are presented in Fig. 7.7. 
These were the first specimens for which the 20:1 lever arm ratio were 
used. Part of the data scatter is related to becoming accustomed to 
using the machine at this ratio. The results from G, H, and I 
laminates are indicated in Figs. 7.8, 7.9, and 7.10. It is not known 
why the data is more consistent for the H specimens than for the 
other two. All three represent quite similar laminates. The data for 
the J specimens is presented in Fig. 7.11 and is seen to be quite 
consistent. At first it was somewhat difficult to obtain good creep 
ruptures for the K specimens. This data is shown in Fig. 7.12. 
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The J and K laminates were composed of only two fiber orientations 
and underwent considerable creep elongation. 

The creep rupture results have been summarized in Table 7.1. 

Based on the best fit lines drawn through the data points, the creep 
rupture strengths for each laminate at 1 minute and at 10,000 minutes 
have been indicated. The percentage decrease in the creep rupture 
strengths over this 10,000 minute time span have also been given. 

While all cases indicate a reduced creep rupture strength at long 
times, the reduction is quite small for several laminates. Interesting- 
v ly the smallest decreases occurred in two laminates with very large 
amounts of scatter, G and I. Strength reductions would possibly have 
been more consistent with the other laminates had less data scatter 
been present. The most important aspect is simply that delayed 
failures have occurred in all the laminates tested. Considering the 
fact that laminates have a design lifetime of many years, the likeli- 
hood for delayed failures for small loads over a long time is quite 
obvious. 

Creep Compliance 

Creep compliance data was taken at several stress levels for most 
laminates. Based on the experimental creep rupture strengths, a 
stress value, was chosen for each laminate which would permit 
compliance testing near the failure strength but would not result in 
failure during the compliance test. Four intermediate, equally 
spaced stress levels, a-| through a^, were also determined. Starting 
with the lowest stress level, short term (36 minute) compliance data 
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TABLE 7.1. Comparison of 1 and 10,000 minute creep rupture strengths 
for the laminates tested. 


LAMINATE 

R(l MIN) 

R(10,000 MIN) 

% DECREASE IN R 

A. 

,C [90/60/ -60/ 90] 2s 

8.86 ksi 

7.94 ksi 

8.5 

D 

[75/45/-75/75] 2s 

9.665 

8.88 

8.1 

E 

[10/55/ -35/1 0] 2$ 

70.35 

67.1 

4.6 

F 

[20/65/ - 25/ 20] 2$ 

57.0 

52.95 

7.1 

G 

[ 90/45/ -45/90] 2s 

20.81 

20.25 

2.7 

H 

[75/30/ -60/75] 2s 

24.4 

22.21 

9.0 

I 

[60/15/-75/60] 2s 

33.03 

31 .98 

3.2 

J 

[15/-75] 4s 

31 .71 

28.97 

8.6 

K 

[30/-60] 4s 

19.45 

16.9 

13.1 
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was taken at 320°F. After an 80 minute or more recovery period, the 
specimen was reloaded at the next highest stress level. At the 
stress level, the test was continued to obtain week long compliance 
data. This testing procedure provided nonlinear behavior information 
based on the short term results, as well as long term compliance data. 
These results have been presented in the subsequent figures, along 
with the compliance predictions for each stress level from the lamina- 
tion program. Because the strain gages on several specimens were 
rendered inoperable before or during the compliance testing, the 
experimental results are incomplete or not present in several figures. 
Compliance testing was begun at a 2 rather than a-j for a few laminates 
because the creep frames could not apply a small enough load. 

Specimens cut at various angles from the [0/30/-30/0] 2s panels 
exhibit much less transient strain than specimens cut from the [0/90]^ 
panel. For example, the 10,000 minute transient strain of the J 
specimen was nearly two orders of magnitude greater than that of the 
E specimen. For specimens containing three or more fiber orientations, 
a vast network of triangular trusses is formed by the fibers. Lamina- 
tion theory assumes that normals remain straight and normal which 
implies that the trusses are effectively "pin connected." These 
triangular structural elements can support the load with the matrix 
material supporting any o 2 or t-j 2 stresses. Compliance predictions 
based on the use of lamination theory will predict an upper bound on 
the compliance for such laminates. This limiting value may be deter- 
mined by allowing the matrix properties to go to zero. 
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On the other hand, laminates with fewer than three fiber orienta- 
tions do not contain the triangular skeleton and do not possess a 
limiting compliance. These laminates rely on the shear panel support 
of the matrix to carry a substantial portion of the load. For these 
reasons ** the time dependent compliance of laminates composed of three 


or more fiber orientations is expected to be much smaller than that of 
laminates composed of less than three fiber directions. It should be 
further noted that for the lamination theory model, an applied creep 
load on the truss-like laminates will result in a "relaxation creep" 
loading for the in-plane matrix stresses within each ply. 

While lamination theory assumes rigid "pin connected" fiber 
trusses, actual laminates undergo interlaminar shear deformations. 
Because these displacements are controlled by the matrix, real speci- 
mens undergo a time dependent relieving of the "pin connection" 
constraint. Thus, experimental compliances are not bounded by the 
compliance asymtote which limits the program predictions. For 
laminates with matrix dominated compliances, the program predictions 
may be quite good. For laminates whose compliance is approaching the 
asymtote, however, the experimental compliances may be significantly 
higher than the predictions. Calculated and actual compliances have 
been schematically illustrated in Fig. 7.13. The use of lamination 
theory places a severe limitation on the prediction of time dependent 
compliances of general laminates. Fiber truss compliances have been 
indicated in the following figures to show how near the predicted 
compliances approach their limiting value. 
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Fig. 7.13 Schematic representation of lamination theory predictions asymptotically approaching the 
limiting compliance value whereas actual compliances are unbounded. 
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Figs. 7.14 and 7.15 show the creep compliance results for the 
C and D laminates. These specimens are from the "old" batch of 
material, while the program predictions were based on the properties 
of the "new" material. Nonetheless, as may be observed, the agreement 
is quite good. The comparisons indicate that the experimental values 
are not as nonlinear in stress as the program predicts. These com- 
pliances are well below the fiber truss limiting value and the pre- 
dicted compliances fan out appropriately . 

The results for laminate E are presented in Fig. 7.16 with the 
fiber truss asymptote indicated near the top of the figure. The 
compliance predictions are seen to be converging as they approach the 
limiting value. While the magnitudes of the compliances agree quite 
well, the experimental short time compliance at the lowest stress level 
is greater than that at the highest stress level. This behavior is 
not consistent with the nonlinear compliance models considered and 
raises some concern. This anomaly was also observed for the F 
specimen orientations, but not in any other laminate. This behavior 
was noted during testing and the test was restarted to verify this 
aspect. Calibration errors associated with changing scales could only 
account for a small amount of such behavior. Our feeling is that this 
phenomena is real and is not simply due to experimental error. Lou and 
Schapery [54] have noted a decreasing compliance after unidirectional 
Gr/Ep specimens have been loaded and unloaded several times. They have 
recommended repeated load application to mechanically condition 
specimens prior to creep compliance testing. Future efforts in this 
area should consider and investigate the possibility of mechanical 
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Fig. 7.14 Comparison of predicted and experimental creep compliance for laminate C ( [75/ 45/-45/75] 2s ) 

at 320°F (160°C) . 
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Fig. 7.15 Comparison of predicted and experimental creep compliance for laminate D ( [75/45/-75/75] 2s ) 
at 320°F (1 60°C) . 
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Comparison of predicted and experimental creep compliance for laminate E ([10/55/-35/10] 2s ) 
at 320°F (1 60°C) . 
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conditioning. Possibly repeated loading at increasingly higher stress 
levels could have some effect on the measured compliances. It is 
felt, however, that the use of successively higher stress levels should 
minimize any repeated loading effects. 

A log-log plot of experimental transient strain vs time is shown 
in Fig. 7.17. The long term data points fall short of the line deter- 
mined from the short term data. Note, however, that these points can 
be made to fall on a straight line by choosing an appropriate value 
of e Q based on the long term data. Figure 7.18 indicates the degree 
of fit for the same data using the Findley approach. e Q appears to 
be a linear relationship, while the values of m are seen to describe 
the characteristic hyperbolic sine function. 

The experimental and predicted compliance results for laminate F 
are shown in Fig. 7.19. Again the anomaly of decreasing compliance 
with increasing stress level should be noted. Because the creep strains 
for both the E and F laminates were quite small, there is some devia- 
tion of the plotted experimental compliance points from a smooth curve. 
Note that the experimental compliance for ag has passed through the 
fiber truss bound on compliance predictions. 

Faulty gages on the G and H laminates prevented obtaining ex- 
perimental compliance data. However, Figs. 7.20 and 7.21 show the 
predicted compliances. As shown in Fig. 7.22, the results for the I 
laminate agree quite well at the low stress levels, although the 
experimental compliance is much larger than predictions at the highest 
stress level. The belief is held that the large measured compliance 
for and a g resulted from matrix cracking in some of the plies. 
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Fig 7 19 Comparison of predicted and experimental creep compliance for laminate F ([20/65/-25/20]2 $ ) 
at 320°F (160°C) . 
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Fig. 7.20 Predicted creep compliance for laminate G [90/45/-45/90] 2s at 320°F (160°C). 
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Fig. 7.21 Predicted creep compliance for laminate H ([75/30/-75/75] 2$ ) at 320°F (160°C). 
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Fig. 7.22 Comparison of predicted and experimental creep compliance for laminate I ( [60/1 5/-75/ 603 2s ) 
at 320°F (1 60°C) . 
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The rapid upturn of the experimental compliance at was due to a 
series of jumps in the strain output. A schematic drawing in Fig. 

7.23 details the nature of these strain steps. The shape of the steps 
tends to indicate that they resulted from localized matrix cracking 
of plies within the region under the strain gages. As some plies 
crack, the remaining laminate became more compliant. The experimental 
transient strains are indicated in Fig. 7.24. 

The compliance results for the J laminate are indicated in Fig. 
7.25. The agreement with the experimental data is considered to be 
very good. The experimental transient strains are shown in Fig. 7.26. 
The results from the application of the Findley procedure to this data 
are shown in Fig. 7.27. The degree of fit is considered to be 
excel! ent. 

Fig. 7.28 presents the K results. Because the J and K laminates 
each contain only two fiber orientations, there is no fiber truss 
limiting compliance. Experimentally, these laminates are character- 
ized by very large transient deformations. Also plotted in this figure 
are the predicted compliances for several stress levels based on 
neglecting the fiber rotations. This laminate underwent the largest 
deformations and fiber rotations. As may be observed, predictions 
are considerably more compliant if significant fiber rotations are not 
accounted for. Unfortunately, the K specimen broke prematurely and 
long term data was not obtained. The transient strain data is pre- 
sented in Fig. 7.29. The Findley results are given in Fig. 7.30 and 
are considered to be very good. 
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Fig. 7.23 Schematic representation of strain jumps for creep loading. 
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Fig. 7.24 Log-log plot of transient strain for Laminate I ([60/1 5/-75/60] 2s ) 
at 320°F (160°C). 
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Fig. 7.25 Comparison of predicted and experimental creep compliance for laminate J ([15/-75] 4s ) 
at 320°F (160°C) . 
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Fig. 7.26 Log-log plot of transient strain for laminate J ([15/75].) at 320°F (160°C) 
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Fig. 7.28 Comparison of predicted and experimental creep compliance for laminate K ( [30/-60] 4s ) at 
320°F (160°C) . 
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Log-log plot of transient strain for laminate K ([30/-60]^ $ ) at 320°F (160°C). 
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Although the Findley procedure was developed for homogeneous, 
isotropic materials, the present endeavor has been to apply this 
approach to the 90° and shear compliances for a unidirectional lamina. 
While this appears to be an appropriate use of the technique, the 
Findley procedure cannot be indiscriminately applied to any arbitrary 
black box material. For example, the predicted compliances of the E 
and F laminates are significantly affected by the limiting fiber truss 
compliance. If the Findley procedure is applied to this model, the 
hyperbolic sine relations are found to be totally inappropriate. 

The implications of applying the Findley approach to the experi- 
mental compliance data of a general laminate have not been fully 
assessed. The procedure, however, was found to be quite appropriate 
in describing the material response, as indicated in the previous 
figures. The power law was also considered to be quite applicable to 
the general laminate response. The results presented herein tend to 
suggest that one could generate carpet plots of the Findley parameters 
for a family of general laminates and such information could be used 
to predict the nonlinear response of arbitrary laminates of the same 
general family. 

Creep Rupture Predictions 

The valid creep rupture data points from Figs. 7.2, 7.5 - 7.12 
are replotted in Figs. 7.31 - 7.39 along with the incremental lamina- 
tion program predictions of the creep rupture strengths. Program 
predictions based on both a = 0.65 and a = 0.80 are presented. 

Program predictions and experimental data are compared with creep 
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Fig. 7.31 Creep rupture data with predictions for laminate C [90/60/-60/90] 2s at 320°F (160°C) 
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Fig. 7.32 Creep rupture data with predictions for laminate D [75/45/-75/75] 2s at 320°F (160°C) 
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Fig. 7.33 Creep rupture data with predictions for laminate E [10/557-35/10^ at 320°F (160°C). 
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Fig. 7.34 Creep rupture data with predictions for laminate F [20/65/-25/20] 2s at 32Q°F (160"C). 
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Fig. 7.35 Creep rupture data with predictions for laminate G [90/45/ -45/90] 2s at 320°F (160°C). 
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Fig. 7.36 Creep rupture data with predictions for laminate H [75/30/-60/75] 2s at 320°F (160°C). 
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Fig. 7.37 Creep rupture data with predictions for laminate I [60/1 5/-7.5/60]2 S at 320°F (160°C). 
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Fig. 7.38 Creep rupture data with predictions for laminate J [15/75]^ at 320°F (]60°C) 
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Fig. 7.39 Creep rupture data with predictions for laminate K [30/-60] 2s at 320°F (160°C). 
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rupture strength predicted by a deformational failure approach which 
assumed that laminate failure would occur when the axial laminate 
strain reached an arbitrarily chosen strain value of 0.5%. Deforma- 
tional failure predictions were based on long term experimental 
compliance data when available. Compliance values predicted by the 
program were used when experimental data was not available. 

Compliance values are a function of stress level because the material 
was nonlinear. The stress values used for the compliances were near 
the creep rupture strengths and are indicated in the figures. 

The results previously given in Fig. 4.5 indicated that a = 0.65 
was an appropriate value for expressing the shear creep rupture 
strength in terms of the 90° strength. In Figs. 7.31 - 7.39, however, 
lamination program predictions based on a = 0.65 fall consistently 
below the experimental data. Also, the results shown in Fig. 4.5, 
indicated that the value of a for the other temperatures tended to be 
greater than 0.65. Based on this evidence and because the program 
predictions all tended to be low, the value a = 0.80 was used to obtain 
another set of program predictions. These values have also been 
indicated in Figs. 7.31 - 7.39. 

The results for the C and D laminates are found in Figs. 7.31 
and 7.32 and were obtained from specimens of the "old" material batch. 
Because input properties used in the program were for the "new" 
material, predictions for the "old" material specimens may not be 
appropriate. Interestingly, predictions based on the two values of a 
converge after a certain length of time and is due to a change in 
the order in which the plies fail. For example, at shorter times in 
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the C specimens, the ± 60° plies are predicted to fail first if 
a = 0.65, but the 90° plies will fail first if a = 0.80. For longer 
times, however, failure is predicted to originate in the 90° plies 
for either value of a, because matrix stresses in the 90° plies relax 
more slowly than in the ± 60° plies. 

The results for E and F laminates are presented in Figs. 7.33 
and 7.34. For the E laminate, a = 0.80 predictions provide a better 
fit to the data than those based on a = .65 but the same is not true 
for the similar F laminate. A possible explanation for this apparent 
discrepancy is that the F specimens encountered much larger bending 
stresses at the grips. Our belief is that end constraint stresses 
have substantially reduced the strength of the F specimens. 

Results for the G, H, and I specimens are given in Figs. 7.35, 
7.36, and 7.37 and, as may be observed, lamination program predictions 
fall consistently below the experimental data. Results for the J and 
K laminates with only two fiber orientations are given in Figs. 7.38 
and 7.39. The a = 0.65 predictions are seen to be low by nearly a 
factor of two, indicating an inability Of the program to predict 
failures in the two fiber orientation laminates. A primary reason 
appears to be based on a great ability of these laminates to continue 
to sustain the creep load after all plies are cracked through the 
thickness. This is examined further in the next section. 

Photographs of Delayed Failure Zones 

The following photographs were taken to illustrate typical 
delayed failure zones in the laminates tested. These plates indicate 
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the manner in which failure occurred and the extent of the damage zone. 
Photomicrographs of typical specimen edges were taken near the 
separation zone and are also presented. The degree of matrix cracking 
in the individual plies should be noted in these plates. Typical C 
and D specimens have been shown in Plates 7.1 and 7.2 respectively. 

Note the very large amount of damage sustained by an E laminate in 
Plate 7.3a. The separation zone for an F specimen is illustrated in 
Plate 7.4a. Note in Plates 7.3b and 7.4b that neither E nor F 
laminates show any apparent matrix cracking away from the actual 
fracture zone. Similar photographs of G, H, and I laminates are pre- 
sented in Plates 7.5, 7.6, and 7.7. The photographs of a J specimen 
in Plate 7.8 are very interesting. Note that all plies were cracked 
through but the laminate was still intact and capable of supporting 
the load. 

The K specimens provided some very interesting details. Plates 
7.9a, 7.9b, and 7.9c illustrate the failure zones for K specimens 
which ruptured at short, medium, and long times, respectively. Note 
the increasing degree of damage with the successively longer rupture 
times. Photomicrographs of these specimens reveals an increasing 
density of matrix cracks with longer rupture times. Several edge views 
of K specimens are illustrated in Plates 7.10 and 7.11. Note the 
more or less regular spacing of the matrix cracks. Stinchcomb, et 
al [71] have discussed such a saturation crack spacing. Although all 
individual plies may be shattered, the laminate remains intact and 
continues to support the load. To ascertain if the cracks actually 
extended across the specimen width, a K specimen was carefully 



Plate 7.1 Typical Creep Rupture Zones for Laminate C [90/60/-60/90] 
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Zone for Laminate E [ 10/ 55 /— 35/ 10 ] 


Plate 7.3b Edge Photomicrograph of Typical E Specimen [ 10/55/-35/10] 
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Plate 7.4a 


Typical Creep Rupture Zone for Laminate F [20/25/-65/20] 2 g • 
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Plate 7.6b Edge Photomicrograph of Typical H Specimen [75/30/-60/75] 2g 
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Plate 7.8b Edge Photomicrograph of Typical J Specimen [15/-75] 




Plate 7.9a Creep Rupture Zone of Laminate K [30/-60] 
t . = .1 min. 


Plate 7.9b Creep Rupture Zone of Laminate K [30/-60], at 18.6 ksi 
t =30 min. S 
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Plate 7.12 Normal View of an Interior -60° Ply of 
a K Specimen [ 30 /— 60 ] Indicating that 
Cracks Extend Across Specimen Width. 




delaminated and a photomicrograph was taken looking down on an 
interior ply. As indicated in Plate 7.12, these cracks do extend 
across the width of the specimen. 

Variation in Ply Stresses 

The variation in ply stresses with time is often thought to be 
quite negligible. The computer predictions of ply stresses, however, 
indicates a considerable decrease in the a ^ and x-^ stresses and a 
corresponding increase in a-j . Figs. 7.40 and 7.41 indicate these 
results for a typical laminate. Fig. 7.40 illustrates the variation 
in a-j , ct^, and x oc ^ for the 90° ply in a G specimen [90/45/ -45/90]2 S 
at 14,500 psi. For this ply, x-^ = 0- Fig. 7.41 illustrates similar 
variation for the 45° ply (or -45° ply). The stresses have been 
normalized with respect to their respective values at t = .01 minute. 
Because the stresses at t = .01 minute were the first iteration in 
the analysis, convergence to correct values may not have been achieved 
This explains why the normalized stress curves do not pass smoothly 
through a value of 1.00 at t = .01 min. The important concept is 
simply that the ply stresses can very significantly due simply to the 
differential creep rates of the various plies. For this laminate, the 
variation exceeded 40% for x-^ of the 45° ply before laminate failure 
occurred. Similar variations were found in the analysis of the other 
laminates. 

Polycarbonate Results 

As mentioned in the previous chapter, some creep yield studies 
were conducted on polycarbonate to investigate certain aspects of 
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Fig. 7.41 Time variation of the normalized ply stresses for 45° ply laminate G 
[90/45/-45/90] Of . at 14,500 psi. * 
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general rate processes. Approximately 80 dogbone specimens were creep 
loaded at 167°F (75°C). Fig. 7.42 illustrates the results for the 
creep tests. The times to yield vary from 1 minute to several days. 

The data is seen to be nearly linear when plotted stress vs log time. 
After this creep yield line for a constant load had been established, 
step-up and step-down tests were conducted to investigate cumula- 
tive damage relationships. These tests were designed so that approxi- 
mately 1/2 of the life to yield was expended at the initial stress 
level prior to the load change. 

Results from the step-down tests show some evidence that the 
linear damage accumulation may be non-conservative for this type loading 
scheme. Yield tended to occur when the linear accumulated yield life- 
time reached about 0.8 rather than 1.0. Of greater interest, however, 
was the evidence that predictions for step-up loadings tended to be 
overly conservative. In fact, it was found that the yield lifetime at 
a given stress level was actually increased if it was first preloaded 
at a lower stress level. Specimens loaded at 5750 psi for 25 minutes 
and then stepped to 6250 psi showed a three-fold increase in lifetime 
at this high stress level over that of a virgin specimen loaded at 
6250 psi. If, however, the specimen was allowed to recover for 25 
minutes following the initial load, the lifetime at the high stress 
level was a fraction of the virgin specimen lifetime. The data taken 
for these step loads was quite minimal and as a result are not pre- 
sented in graphical form. More extensive testing should be conducted 
to verify these indications. 




Fig. 7.42. Creep yield of 
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The conclusions, however, is that the linear damage accumulation 
rule is not always appropriate for polycarbonate. The loading 
sequence can be very important because, apparently, some type of 
beneficial aging process takes place at the low stress level. Because 
this improvement is negated when the specimen is allowed to recover 
prior to applying the high load, this beneficial aging must be a 
mechanical or stress effect rather than of thermal origin. 

Physical Aging Effects in Graphite/Epoxy 

Mechanical aging effects similar to those in the polycarbonate 
have also been observed in the current testing of Gr/Ep specimens. 

Again it was found that a creep preload tended to increase the ultimate 
breaking strength of a specimen. If a specimen had not broken under 
creep loading after about a week, the load was gradually increased by 
the addition of one pound weights to the weight pan until the specimen 
broke. This breaking strength was considerably higher than the static 
strength of a virgin specimen. 

The results of the breaking strengths of the specimens for the 
preloaded specimens are given in Table 7.2. The creep load and its 
duration are given along with the breaking strength. These strength 
values are compared with the estimated creep rupture strength of the 
given laminate at 1 minute. These values were obtained from the best 
fit lines in the creep rupture figures. As may be seen there was 
about a 15% increase in static strength of the specimen if it had been 
preloaded. The K specimens, composed of only two fiber orientations, 
showed an even greater increase in strength and underwent the largest 
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TABLE 7.2. 

Increase 

preload. 

in static 

strength due to 

mechanical 

aging of a 

SPECIMEN 

PRELOAD 

(KSI) 

TIME AT 
PRELOAD 
(MIN) 

BREAKING 

STRENGTH 

(KSI) 

R (1 MIN) 
(KSI) 

% INCREASE 

E-14 

67.0 

9850 

79.0 

70.4 

12.2 

F-2 

44.0 

11440 

68.0 

57.0 

19.3 

F-n 

54.0 

10000 

64.8 

57.0 

13.7 

F- 

54.5 

9400 

70.8 

57.0 

24.2 

1-8 

31 .0 

9870 

38.2 

33.0 

15.8 

J-2 

25.7 

3100 

35.9 

31 .7 

13.2 

J-3 

27.0 

2640 

36.2 

31 .7 

14.2 

J-4 

28.2 

2570 

38.6 

31 .7 

21 .8 

J-5 

29.0 

9837 

37.4 

31 .7 

18.0 

J-19 

29.2 

9533 

35.2 

31 .7 

11.0 

K-6 

18.0 

11360 

26.4 

19.4 

36.1 

K-7 

18.8 

9790 

25.9 

19.4 

33.5 



198 


deformations. This strengthening effect may have been due in part to 
the smaller fiber angles which result from the scissoring action 
associated with large deformations. 

Strengths tended to approach that of a virgin specimen if the 
specimens were allowed to recover several days after the initial creep 
loading. 

Despite the major differences between the ductile polycarbonate 
and brittle epoxy matrix, significant similarities exist between these 
materials. In particular, the mechanical aging phenomenon in both 
materials seem to be closely related. Both materials appear to exhibit 
a reversible strength improvement from a moderate level creep load. 
Recovery tends to return the material back to a quasi -virgin state. 

Accuracy of Predictions 

In general the predicted compliance values are considered to be 
quite promising when compared with the experimental data. The short 
time compliance predictions are all quite accurate. The compliance 
predictions at long times are not as correct, perhaps because of errors 
in the power law model. Except for the erratic results of the E and 
F laminates, the agreement of the transient nonlinear effect is fair. 
The predictions for laminate compliances approaching the fiber truss 
asymptote are significantly in error. This behavior is a characteris- 
tic of the lamination theory model used herein. If the approach could 
be modified to relax the interlaminar shear deformation constraint, 
better compliance predictions might be obtained without resorting to 
a more costly finite element analysis procedure. The authors are 
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encouraged by the accuracy of the results. 

The creep rupture strength predictions based on apparent 
strength properties of the unidirectional material are low for all 
laminates. Although conservative predictions are good, the current 
results are often overly conservative. The failure model is believed 
to be responsible for these discrepancies because of the assumption 
that a ply ceases to support any o ^ or stresses once ply matrix 
failure has occurred. As a result, the laminate loses all predicted 
load carrying capabilities when the matrix in all plies has broken. 

The experimental results indicate that this assumption may not be 
appropriate. 

Often the failure of several plies within a general laminate is 
not felt to necessarily result in total laminate failure. In the 
current data there are some indications that this is correct. The 
load-deflection curves for several laminates tested at a constant 
crosshead speed exhibited plateaus in the load response such as the 
results for the D laminate which have been illustrated in Fig. 7.43. 
Jones [48] has indicated such plateaus may represent the failure of 
certain plies within the laminate and he has idealized such behavior 
with a laminate model composed of springs in parallel. If each spring 
represents a particular ply, failure of a certain spring would result 
in a horizontal shift in the load-deflection curve. The results given 
in Figs. 7.22 and 7.23 for specimen I indicate probable matrix cracking 
at the stress level, yet the specimen supported an even higher 
stress at a subsequent load level. The edge view photomicrographs 
indicated that several laminates were significantly cracked even at 
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Fig. 7.43 Load-deflection curve for a D specimen ([75/45/-75/75] ? ) 
in constant crosshead loading. 
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large distances from the failure zone. Thus there was substantial 
evidence that first ply matrix cracking does not necessarily result 
in immediate failure of the total laminate. 

The program predictions for the laminates tested, however, indi- 
cated that matrix failure in the plies of a single orientation always 
resulted in almost immediate global laminate failure following the 
first ply failure. These predictions were based on eliminating all 
(?£ and stiffness contributions from plies in which matrix cracking 
had been predicted. Based on the experimental evidence, more accuracy 
may have been obtained by allowing broken plies to remain partially 
effective. Determination of appropriate methods to achieve this 
effect, however, could be quite difficult. 

The assumption that a laminate could not support any load if all 
plies had cracked was also not borne out experimentally . Indeed, the 
individual plies in the J and K laminates are extensively shattered, 
and yet the laminate continued to support the applied creep load. 

The failure model does not account for the strength contributions made 
by the intact regions between the matrix cracks within the plies. 
Interlaminar shear stresses can continue to support an applied load 
via a contorted load path. Unless this strength contribution can be 
accounted for, failure predictions should be expected to be low. 

A crack may easily propagate over the cross-sectional area and 
result in a clean break in unidirectional laminates such as those used 
to determine the material failure properties. However, a flaw cannot 
propagate as a smooth break for general laminates. The actual damage 
region for a general laminate may be quite extensive, as indicated in 
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the photographs of the laminate failure zones. Flaws must also grow 
between the plies to provide the delamination necessary for separation 
contributing to the insensitivity of laminated fiber reinforced 
materials to internal flaws. There may be inherent difficulties in 
applying the flaw sensitive unidirectional failure properties to 
general laminates which are not as susceptible to global failure 
resulting from a flaw in a single ply. 

For these reasons, the laminate failure model used is felt to 
be overly conservative. Two other aspects tend to make our model non- 
conservative, but apparently do not completely compensate for the 
problem. Lamination theory ignores all interlaminar stresses which 
are known to exist at the free edges and constraint stresses intro- 
duced by the grips have also been neglected. These two factors produce 
stress states that are greater than those accounted for by the lamina- 
tion theory model. Development of a more accurate failure model to 
account for all the above effects could be quite difficult. 

A prominent discrepancy exists for all laminates between the rate 
at which the creep rupture strength predictions decrease with time and 
that indicated by actual creep rupture data. This difference cannot 
be explained by lamination theory shortcomings. In fact, lamination 
theory predicts a greater decrease in matrix stresses for laminates 
near the fiber truss compliance limit than actually occurs. This 
should have the effect of decreasing the predicted rate of creep rupture 
strength reduction. 

Failure properties determined from non-post-cured unidirectional 
creep rupture data may have contributed' to the rapid decrease of 
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predicted laminate creep rupture strengths. As obtained from Eqns. 

4.6b and 4.6c, the decrease of Y(t ) and S(t r ) from 1 to 10,000 
minutes is 32% as opposed to the 3 to 13% decrease obtained for the 
general laminates tested, as indicated in Table 7.1. Obviously, the 
use of Eqns. 4.6b and 4.6c should result in rapidly decreasing strength 
predictions. Recently Eqn. 4.6b was noticed to have been based on 
Griffith's [38] best fit line for the 90° specimens at 320° F. This 
line does not agree with the best fit line drawn in Fig. 4.2. Our 
best fit line indicates a 26% strength reduction rather than the 32% 
for the 10,000 minute time span. Use of this corrected best fit line 
would result in a small improvement over the results obtained using 
Eqn. 4.6b. 

Again, difficulties are apparent for determining functional 
relationships from a minimal amount of widely scattered data. If the 
indicated strength reductions for the unidirectional and general 
laminates are correct, the rate discrepancy between experimental and 
predicted strengths may allude to errors in the cumulative damage law 
or the constitutive model. Our lamination program predicts that the 
ply stresses decrease with time. The polycarbonate creep yield data 
indicated that the linear cumulative damage law was non-conservative 
for step-down loading. Although not verified for Gr/Ep, the implica- 
tion of the latter is that the predicted strengths would decrease more 
rapidly than the experimental values which is in general agreement 
with the comparisons previously given in Figs. 7.31 - 7.39. 

One other possible explanation for the rate discrepancy is that 
the unidirectional creep rupture data obtained by Griffith [41] was 
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for specimens which had not been postcured. Creep ruptures for 
several postcured 60° specimens were obtained during the current work 
and have also been illustrated in Fig. 4.2. The best fit line through 
this postcured data indicates only a 9% decrease in the 60° creep 
rupture strength over the 10,000 minute time span rather than the 26% 
decrease for non-postcured specimens. The implication is that the 
rates of strength reduction may be considerably smaller than for non- 
postcured specimens. The data scatter prevents substantiation of this 
claim. Because the failure model was based on no postcure yet the 
general laminates tested were postcured, this is a possible explanation 
for the prediction discrepancies. 

Interestingly Figs. 7.31 - 7.39 indicate that the curves based 
on failure at .5% strain are often in better agreement with the 
creep rupture data than the program predictions. In particular, the 
slope of the predicted and experimental creep rupture strengths are 
often quite similar. This was true for all cases except the two- 
fiber orientation laminates, J and K. One shortcoming of such a 
procedure is that the creep rupture strength is predicted to decrease 
more rapidly at larger values of (log t r ). Physically, however, one 
would expect creep rupture curves to flatten out in order that a zero 
stress would not result in a finite time to rupture. 

Grip Constraint Stresses 

Because of the shear coupling effect in unbalanced laminates, an 
axial stress results in shear as well as normal deformations. A 
complex shear and bending moment stress state is developed in such 
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specimens because of the constraints imposed by fixed grips. A static 
finite element model was used to determine the variation in stress 
state throughout each laminate but difficulties were encountered in 
trying to model the boundary conditions. The procedure was abandoned 
in favor of a closed form solution provided by Pagano and Halpin [56] 
and based on some simplifying assumptions about the boundary condi- 
tions. Their equations were used to predict the axial stresses at the 
grip on each edge of the specimen. These locations represent the 
maximum deviations of axial stresses from the nominal value. Nor- 
malized results have been given in Table 7.3. A general correlation 
between large deviations from the nominal stress value and the likeli- 
hood of failures near the grip has been noted. Pagano and Halpin' 
have also indicated that the apparent stiffness measurements may be 
different than the actual stiffness because of the shear coupling. 

The apparent modulus is given in terms of the actual modulus and a 

parameter z 
* 

E = E ? 

XX XX 

Values of t have also been tabulated in Table 7.3. The parameter z is 
highly dependent on the specimen aspect ratio. For the specimen length 
used, the error in the apparent modulus is quite small except for the 
J laminate. The values in Table 7.3 are based on no lateral motion 
of the grips. Because of the small lateral constraint on the creep 
frames, the actual error for our work is considered to be even less 
than indicated in the table. 



TABLE 7.3. Effect of grip constraint on laminate 
stresses and apparent modulus. 


SPECIMEN 

a /a 

max' nom 

a • /a 
min nom 

4 

D 

.9034 

1 .091 

1 .0015 

E 

.9774 

1 .022 

1 .00009 

F 

.935 

1 .062 

1 .00069 

H 

.9217 

1 .074 

1.00098 

I 

.8539 

1 .133 

1 .0033 

J 

.6088 

1 .306 

1 .0219 

K 

.7986 

1 .177 

1 .0061 



Chapter 8 

SUMMARY AND RECOMMENDATIONS 

An incremental numerical scheme based on lamination theory was 
developed to predict the time dependent response of general laminates 
composed of orthotropic laminae. The procedure uses a nonlinear 
compliance model based on an extension of a technique due to Findley 
to a biaxial stress state in a fiber reinforced material. The octa- 
hedral shear stress, based on the stress state in the matrix, was used 
as the nonlinearizing parameter in order to provide an interaction 
effect among the stresses. The time variation of the compliances was 
assumed to obey the power law for creep. Findley's modified super- 
position principle was used to determine the strains resulting from a 
time varying stress state. Also incorporated into the numerical pro- 
cedure was a plywise failure model based on a modification of the 
Tsa i -Hi 1 1 criteria which accounts for time dependent creep rupture 
strengths. Thus, the 90° strength rupture time was assumed to increase 
exponentially with decreasing applied stress. The shear strength as 
required for the Tsai -Hi 11 approach was assumed to be a constant 
fraction of the 90° creep rupture strength. A linear damage accumula- 
tion rule was used to determine the life expenditure at each time step. 
The procedure incremented through time until all plies were deter- 
mined to have failed. The material properties used in the analysis 
were determined from tests on unidirectional specimens. 

Delayed failures were produced in a variety of different lami- 
nates. The specimens ranged from matrix dominated layups to those in 
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which the load axis was just slightly off-axis from the predominant 
fiber direction. Both balanced and unbalanced laminates were tested. 
Creep rupture data indicated that at 320°F stress levels 3 to 13% 
below the short time strengths could result in delayed failures within 
one week. Considerable scatter characterized the creep rupture data 
for several laminates studied. One of the panels used contained plies 
with crooked fibers and these were believed to have precipitated pre- 
mature failures. Some specimens exhibited considerable thickness 
variations and presented difficulty in determining effective cross- 
sectional areas. Furthermore, particularly in the unbalanced laminates, 
the stress state varied considerably along the specimen length. These 
factors could have led to even more significant scatter than that 
normally observed for creep rupture of homogeneous, isotropic 
materials. 

Based on the failure model developed from uniaxial tests, 
computer predictions of creep rupture strengths were presented for 
each laminate tested. In general, the computer strength predictions 
were lower than the experimental creep rupture data. A possible reason 
may be due to the assumption that total laminate failure occurred when 
the matrix in each ply failed. Photomicrographs of the specimen edges 
indicated that plies could continue to be partially effective although 
the matrix of each ply was cracked through. Laminates could remain 
intact and thus continue to support the applied load even after the 
matrix of each ply had cracked. The conservativeness of this assump- 
tion is partially offset by lamination theory's neglect of the effects 
of interlaminar stresses. Also, the shear and bending stresses induced 



by the grip constraints in unbalanced laminates were not considered. 


The predicted creep compliance strengths decreased more rapidly 
with time than was indicated by the experimental data. A possible 
explanation for this discrepancy is that the failure model was based 
on creep rupture data of non-postcured unidirectional specimens, 
whereas the laminates tested were all postcured. Several postcured 
60 ° specimens were tested and the results indicated that the postcure 
process did not change the basic strength substantially. However, 
postcured specimens did tend to exhibit a smaller decrease in creep 
rupture strength with time than those that had not been postcured. 

Creep compliance data was taken for several laminates and was 
compared with predictions made by the lamination program. The assump- 
tion that no interlaminar deformations exist, results in an upper 
bound on the predicted compliances for laminates composed of three or 
more fiber orientations. The fibers in these laminates produce a 
triangular truss network which fixes the maximum predicted laminate 
deformation unless interlaminar displacements are permitted. For 
laminates which did not possess a limiting compliance value, or for 
those in which the measured compliances were well below this 
asymptotic value, the predicted compliances agreed quite well with 
the experimental values. The upper bound provided by the fiber- truss 
network, however, proved to be a severe shortcoming of the applica- 
tion of lamination theory to certain laminates. A finite element 
approach could have avoided this limitation, but, as discussed earlier, 
would likely have been very costly to implement. If the constraint 
caused by the assumption that no interlaminar deformation exists could 
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be partially relaxed, the lamination theory approach might be more 
successful . 

A mechanical aging phenomena has been observed in both Gr/Ep 
and polycarbonate. Similarities of the rate processes in these two 
very different materials have been noted, particularly with regard to 
the mechanical aging phenomenon. 

Our investigation has been based on integrating a variety of 
appropriate concepts into a single, procedure to predict the compliance 
and delayed failures of general laminated composites composed of a 
fiber reinforced nonlinear viscoelastic matrix. In developing the 
procedure, aspects have been drawn from a variety of seemingly un- 
related areas such as: lamination theory, orthotropic viscoelasticity, 

static laminate failure, creep of metals, cumulative damage laws, 
metal plasticity, numerical procedures, and previous experimental 
data for the unidirectional material. As with any predictive tech- 
nique, the overall accuracy can be no better than that of its 
component assumptions. While the predictions have borne out certain 
trends reasonably well, they have not been acceptable in several ways. 
In pursuing any future analytical treatment, it is proposed that many 
of the assumptions used within this work should be verified or modi- 
fied as more data is obtained. 

Experimental Recommendations 

The unidirectional properties used in the current analysis were 
determined from existing data, which was at times widely scattered 
and even somewhat inconsistent. Most of these results were obtained 
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prior to the realization that there was a significant postcuring 
effect. As such, much of the existing data represents a rather 
haphazard degree of postcuring which of necessity took place during 
each creep test at elevated temperature. Accurate expressions for 
material response are essential for developing predictive techniques. 

Many aspects concerning the characterization of Gr/Ep have been 
revealed by the current project during the past several years. With 
these new concepts in mind, perhaps the time has come to systematically 
repeat much of the testing. A sufficient quantity of material should 
be obtained so that all specimens for the foreseeable testing program 
could be cut from the same batch of material. Care should be taken 
that all specimens receive the same postcure treatment and contain 
similar moisture contents. The recommendation is made that mechanical 
conditioning effect on compliance, as discussed by Lou and Schapery 
[54], be investigated. Testing should be tailored to provide informa- 
tion required to characterize the most appropriate models for 
characterizing the various facets of material behavior which contribute 
to the creep rupture process. 

The use of compliance and creep rupture data obtained at several 
temperatures to predict material response at other temperatures has 
been advanced in previous work [38,82]. The present work has been 
conducted primarily at only one temperature. The feeling is held that 
in developing compliance, creep rupture, and numerical models for 
the Gr/Ep material studied herein, concentration on a single temperature 
is the best procedure. Once acceptable techniques have been developed 
at a particular temperature, generalizations can be made for the 
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thermal dependence. 

Recommendations for the Compliance Model 

The Findley approach to nonlinear compliance has been found 
to be quite appropriate in matching data from unidirectional and 
general laminates. It is believed that the power law provides an 
adequate description of the compliance time dependence. Determining 
the exponent from short term data is not generally valid for pre- 
dicting long term compliance. The short term data generally dictates 
a higher value of the exponent than is indicated by long term response. 
While this points to the inappropriateness of the power law, it is 
believed that an exponent obtained from long term data can give 
quite accurate predictions over a very large time span. The power 
law is not exact, but it is felt to be the best simple approach 
available at this time. A further problem associated with using 
short term data to evaluate the exponent is the errors introduced by 
the singularity problem. There is some evidence that the power law 
exponent may vary with the applied stress level. It is recommended 
that week long compliance data be taken for several stress levels on 
10° and 90° specimens. This long term data should provide a good base 
to establish the dependence of the power law exponent on the applied 
stress level. 

It has been assumed that the averaged octahedral shear stress 
in the matrix serves as the nonlinearizing effect for the compliance. 
There is no verification for this assumption on the Gr/Ep system 
studied. More extensive compliance testing at several off-axis angles 
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could provide information on this or other interaction terms. 

The exponents for Griffith’s [38] 30° and 60° long term com- 
pliance were significantly higher than those for the 10° and 90° 
tests, as indicated in Chapter 7. The transformation equations cannot 
predict these differences. It is recommended that these tests be 
repeated to determine the accuracy of the results. If the previous 
data is correct, a substantial rework of the compliance model and/or 
transformation equations would be required. The proposed tests could 
be the same as indicated in the previous paragraph. 

The modified superposition principle used to account for a 
varying stress level was felt to be adequate for the small stress 
variations associated with creep loaded laminates. However, this 
approach is not believed to be accurate for a more general stress 
variation. Because of its more general nature, the Schapery pro- 
cedure should provide a better treatment of this aspect. While this 
technique does require more extensive testing, it could provide a more 
unified compliance model. The hyperbolic sine terms of the Findley 
procedure have been found to provide good agreement with the experi- 
mental data. Perhaps these terms could be incorporated into the 
Schapery procedure to assist determination of the stress dependent 
functions. It has been proposed that multiple step loads be applied 
to 10° and 90° specimens to determine the appropriateness of the 
modified superposition procedure used herein. Such information could 
also prove useful in evaluating the performance of the Schapery 
procedure for our material system. 
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Recommendations for the Failure Model 

Evaluation of meaningful functional relationships for creep 
rupture of unidirectional specimens is hampered by the severe 
data scatter. Accurate property descriptions are imperative when 
these values will be used to predict failure in other configurations. 

It is proposed that an extensive testing program be conducted to obtain 
creep rupture data for the unidirectional material at several off-axis 
angles. This could establish the appropriateness of the modified 
Tsai-Hi 1 1 failure criteria used herein or could lead to the develop- 
ment of a better creep rupture criteria for unidirectional specimens. 

The cumulative damage law used is believed to be inaccurate. 
Cumulative damage in metals is currently supporting a great deal of 
investigation. Perhaps this technology could provide better analysis 
procedures. The indications of a mechanical aging phenomena should 
be investigated as their effects may continue to frustrate applications 
of cumulative damage laws. 

It is possible that other time dependent failure models would 
prove to be more accurate. Three basic approaches have been widely 
used in the past to predict delayed yield and rupture. By far the 
most attention has been given to the stress approach as used herein. 

The deformation and energy approaches have also been studied. These 
methods may prove more adaptable to the damage accumulation concept. 

Brinson [9] has advocated the use of a deformational failure 
criterion in a ply by ply fashion to predict delayed failures of 
general laminates. Specifically, he has postulated that for a uni- 
directional material, the shape of the creep rupture curve is the same 
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as that of the inverse of the creep compliance. Short term compliance 
data at several temperatures could be shifted to obtain a master 
curve of the creep compliance. The implication of this approach is 
that this compliance function could be inverted and normalized to 
obtain a failure master curve. If this procedure could be verified, 
a substantial reduction in required creep rupture testing would result. 
A simple extension of this procedure to a general biaxial stress state 
in a lamina is not apparent, however. 

Bruller [10,11,12] has advocated the use of an "energetical 
limit" in predicting the viscoelastic yield of polymers. His approach 
is based on the Reimer-Weissenberg Theory which states that visco- 
elastic yield will occur when the stored deviatoric energy reaches a 
critical value known as the resilience - a material property. While 
this approach seems quite accurate for predicting yield in polymers, 
it is limited to linear viscoelastic materials which can be modeled by 
a generalized Kelvin element. 

One disturbing feature about the approach used herein is that 
the compl iance and failure phenomena are addressed by unrelated models. 
It would seem advantageous to provide a common basis for these two 
aspects. Perhaps the use of a deformation failure criteria could 
unify the treatment. It should be noted that the creep rupture theory 
indicated by Eqn. 4.2 is a power law similar to that used for the 
compliance model. With these approaches, however, one has still not 
properly addressed the superposition effect or cumulative damage of 
a time-varying stress state. It was originally thought that the 
Tobol sky-Eyring Reaction Rate Equation could be adapted to treat this 
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problem. This equation relates the rate at which atomic bonds are 
broken, and is given by 


dN 

dt 


kT 


exp 


F 

RT 


2 sinh 


fwj 

RT 


where 

N = number of bonds/unit cross section 
&F = free energy of activation 
T = temperature 
k = Boltzmann's constant 
h = Planck's constant 
R = gas constant 

W = work done on a single bond applied by a. 

It has been widely studied by a number of investigators [1,2,3, 
42,46]. Because the number of remaining bonds represents a state 
variable, it is possible that this could provide a unified approach to 
compliance and failure under a biaxial, time-varying stress state. 

Recommendations for the Numerical Procedure 

The most significant flaw in the compliance predictions is the 
upper bound imposed by the fiber truss work. This effect results 
from the manner in which the individual ply stiffnesses are combined 
according to lamination theory. By artificially reducing the lateral 
stiffness, one could relax the constraint imposed by the assumption 
that normals remain straight and normal. If this does not prove 
feasible, one may be forced to use a finite element approach. 
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The computer predictions for strength are low. This could 
possibly be remedied by allowing "failed" plies to remain partially 
effective. As illustrated in the damage zone photographs, the 
failure mechanism for general laminates is quite complex. Accounting 
for these intricate processes will be very difficult. Some of the 
procedures used to predict static failures of general laminates are 
not applicable to delayed failures. Considerable work remains to 
be done in this area. 

It is expected that a finite element approach will eventually 
be desired to analyze general laminates. More valid models for 
material response will warrant the improved accuracy afforded by the 
finite element method. A three-dimensional approach could model the 
effects of the end constraints as well as the interlaminar stresses. 

It might be possible to develop a two-dimensional element to represent 
the specimen cross-section. This could account for the interlaminar 
stresses without requiring so much storage. 

Conclusion 

The procedure developed is not the final answer to analyzing 
time dependent behavior in laminated composites. Nonetheless, it is 
felt that this work can help pave the way for developing better tools 
towards this end. The author envisions developing better compliance 
and failure models which could eventually be incorporated into a more 
accurate finite element solution scheme. 
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Appendix A 

LEAST SQUARES HYPERBOLIC SINE FIT 

The following procedure was used to obtain a least squares fit 
of a hyperbolic sine function 

y = a sinh(x/b) 

to a given set of data (x^, y. ) i = l,2,...,n 

Let 

e i = y - y i = a sinh(x/b) - y. 
and 

n p o o p 

I = z e. = a z sinh (x-/b) - 2a z sinh(x./b)y. + z y. 
i=l 1 1 iii 

Now let 

O 

a = z sinh(x.j/.b)y ./£ sinh (x^./b) 
and 

F(b) = [z sinh(x i /b)y i ][z sinh(x.j/b) cosh(x./b) (x^b 2 )] 

- [z sinh 2 (x i /b)][z coshfxj/b) (x i /b 2 )y i ] = 0 

Roots to this equation will render I stationary with respect to a and 
b. Solving for the roots is not as easy task because of the shape of 
the F function. The secant method (modification of the Newton - Raphson 
procedure which uses derivatives approximated by successive secants) 
was tried, but obtaining a starting value which did not result in an 
overflow or diverge to infinity was quite difficult. Use of the 
bisection method was found to yield an appropriate initial guess. The 
combination of these two techniques provided an adequate solution scheme. 
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Appendix B 

EFFECTS OF LEVER ARM OSCILLATION AND ROTATION 

In creep rupture testing, the specimen may linger on the verge 
of failure for some time prior to actually rupturing. When super- 
imposed on the applied load, small additional stresses such as machine 
vibration may be sufficient to trigger a premature failure. Such 
oscillations should be minimized to insure accurate results. 

Load train oscillation is an annoying feature of the ATS Lever 
Arm Test Frame at the 20:1 arm ratio. The specimen and load train 
mass act as a simple spring and mass system to result in the 
vibration problem. Aside from the friction of the knife edge 
supports, the only significant attenuation is the internal damping of 
the specimen. Applying the load by lowering the weight elevator 
resulted in large oscillations. The load may be applied more gradually 
by using the variable speed drive for the crosshead while set at a slow 
rate. This significantly reduces the system oscillation, but does 
not eliminate them. Measurements were taken to evaluate the dynamic 
effects. At the 20:1 ratio, a pan load of 150 lbs resulted in a 
stress of 70,500 psi on specimen E-l . An effective weight for the 
entire load train was taken to be 165 lbs. A dial indicator mounted 
on the test frame was used to monitor the weight pan deflection. 

By using a gradual application of the load as described above, the 
maximum load displacement did not exceed 5 = .01". The oscillation 
frequency was about 2 Hz. This resulted in a sinusoidal dynamic load 
of amplitude .34 lbs or a load variation of .22%. Because the stress 
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levels for our creep ruptures were all within several percents of one 
another, this may not be a negligible quantity. Care should be 
taken to minimize any machine vibration. 

Occasionally a concern is raised regarding the variation in 
applied stress of the lever arm creep machines as the specimen 
elongates and permits the lever arm to droop. Measurements taken on 
the Budd machine were used to determine the maximum error associated 
with this aspect. For a specimen elongation of 0.5", the lever arm 
rotates 10° and the load drops 5". Based on simple geometric con- 
siderations, the resulting deviation in the applied stress is only 
0.03%--a negligible amount. As the effective lever arm of the applied 
load decreases with arm rotation, the effective lever arm of the 
specimen load train decreases by an almost identical proportion. 
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